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Abstract 


A  neutron  time-of-f light  system  using  a  5.5  MeV  Van  de  Graaff  with 
a  Mobley  compression  magnet  was  set  between  the  target  and  liquid  helium 
polarimeter  for  the  measurement  of  fast  neutron  polarization.  The  polari¬ 
zation  system  was  calibrated  with  neutrons  from  the  D(d,n) 3He  reaction  at 
deuteron  energies  of  3.0,  4.0,  5.0  and  6.0  MeV  and  at  a  neutron  emmision 
angle  of  45°(c.m).  Good  agreement  was  found  between  the  present  measure¬ 
ments  and  those  of  Dubbeldam  and  Walter. 

The  polarization  angular  distributions  of  neutrons  to  the  ground, 

1.38  and  1.96  MeV  states  of  29P  from  the  reaction  28Si(d,n)29P  were 
measured  at  a  deuteron  energy  of  5.0  MeV.  The  polarization  angular  dis¬ 
tribution  corresponding  to  the  ground  state  (£p=0)  group  is  in  good  quali¬ 
tative  agreement  with  the  derivative  rule.  Spin  dependent  effects  are 
observed  in  the  polarization  distributions  of  the  1.38  (j=£  -  -j)  and  the 
1.96  (j=£p+  y)  MeV  states. 

The  DWBA  code  written  by  B.E.F.  Macefield  was  used  in  the  calculation 
of  neutron  polarization.  Fair  agreement  between  the  DWBA  prediction  and 
the  polarization  measurements  to  the  ground  state  of  28P  is  found,  but 
there  is  very  poor  agreement  for  the  1.38  MeV  state.  As  for  the  1.96  MeV 
state  there  is  some  agreement  in  the  shape  of  the  polarization  distribution 


in  the  forward  angles. 
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INTRODUCTION 


In  the  study  of  neutron  producing  reactions,  although  both  the  incident 
beam  of  particles  and  the  target  are  unpolarized,  the  outgoing  neutrons 
are  in  general  found  to  be  polarized.  This  is  due  to  the  spin-orbit  coup¬ 
ling  in  nuclear  reactions.  Since  the  polarization  vector  of  the  outgoing 
neutrons  is  normal  to  the  reaction  plane  and  the  neutron  has  an  intrinsic 
spin  of  1/2,  the  neutron  spin  orientations  can  only  be  up  or  down.  If  in 
a  beam  of  neutrons,  the  number  with  spins  up  is  N+  and  those  with  spins 
down  is  N_,  the  polarization  of  the  beam  can  be  defined  as: 

N+  -  N_ 

P  =  - 

N_j_  +  N_ 

If  N+  is  equal  to  N_,  the  polarization  is  zero.  If  either  N_  or  N+  is 
zero,  the  polarization  is  +  1. 

Schwinger  (Sc  46,  Sc  48)  was  the  first  to  indicate  the  importance  of 
polarization  studies  in  nuclear  research.  Lepore  (Le  50)  made  detailed 
calculations  based  on  Schwinger's  suggestions  and  Wolfenstein  (Wo  49,  Wo  54, 
Wo  56)  made  extensive  studies  in  the  case  of  nucleon-nucleon  scattering. 
Blin-Stoyle  (B1  51)  was  the  first  to  make  the  observation  that  polarized 
neutrons  should  be  a  frequent  product  of  nuclear  reactions.  Blatt  and 
Biedenharn  (B1  52)  were  first  to  give  a  complete  account  of  the  polarization 


angular  distribution. 


- 

. 
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The  picture  of  particle  polarization  from  a  stripping  reaction  was 
first  put  forward  by  Newns  (Ne  53)  in  terms  of  nuclear  absoptive  distortion. 
Since  the  Butler  plane  wave  theory  (Bu  51)  does  not  account  for  polarization, 
the  distorted  wave  Born  approximation  (Au  63)  is  introduced.  Neglecting 
the  spin-orbit  term  in  the  optical  potentials,  the  upper  limit  for  polari¬ 
zation  is  33-1/3  %  and  for  the  case  where  the  orbital  angular  momentum 
transfer  of  the  captured  particle  is  zero,  the  prediction  for  polarization 
is  zero  (Hu  58,  Bi  60).  Since  polarization  measurements  have  shown  these 
to  be  not  true,  spin-orbit  potentials  are  introduced  in  the  distorted  wave 
Born  approximation.  It  has  been  shown  that  to  the  first  order  the  angular 
distribution  is  independent  of,  while  polarization  depends  linearly  on,  the 
spin-orbit  potentials  (Jo  62),  so  polarization  measurements  should  be  a  more 
sensitive  way  of  testing  the  strength  of  spin-orbit  potentials  in  a  nuclear 
reaction.  With  the  availability  of  high  speed  and  large  memory  computers, 
programs  (Ba  61,  Sa  64,  Ma  64,  To  66)  have  been  written  for  distorted  wave 
Born  approximation  calculations  with  spin-orbit  potentials. 

The  object  of  this  thesis  is  to  report  on  the  construction  and  applica¬ 
tion  of  a  liquid  helium  polarimeter  and  the  comparison  of  neutron  polari¬ 
zation  measurements  from  the  reaction  28Si(d,n)28P  with  the  distorted  wave 
Born  approximation  predictions. 

There  has  been  a  conflict  in  sign  convention  in  discussions  of  polari¬ 
zation.  The  Basel  sign  convention  (Ba  60A)  is  adopted  in  this  thesis.  This 

->■ 

prescribes  a  positive  sign  of  polarization  in  the  direction  kin  x  kout, 

-4- 

where  k^n  and  kout  are  the  momenta  of  the  incident  beam  of  particles 
and  the  outgoing  particles,  respectively. 
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CHAPTER  I 

LIQUID  HELIUM  POLARIMETER 

1-1  Introduction 

Since  neutrons  produced  from  a  nuclear  reaction  are  usually  polarized 
and  conventional  detectors  cannot  distinguish  the  spin  orientations  of  the 
neutrons  detected,  special  arrangements  have  to  be  used  to  measure  polari¬ 
zation.  While  experiments  using  a  polarized  beam  of  particles  or  a 
polarized  target  have  been  performed,  elastic  scattering  from  nuclei  is 
most  commonly  used  in  the  study  of  polarization.  The  various  techniques 
in  the  measurement  of  polarization  have  been  described  by  Haeberli  in 
Chapter  V.G.  of  Fast  Neutron  Physics  Part  II  (Ha  63)  and  the  following 
description  will  be  devoted  entirely  to  the  system  adopted  at  this  laboratory. 

Fig.  1-1  shows  a  typical  arrangement  in  a  polarization  experiment.  If 
an  incident  beam  of  particles  with  momentum  kj  hits  a  target,  the  polarization, 
Pi,  of  the  neutron  emitted  at  an  angle  6i  is  given  by: 

P1  =  P1  nl 

where  ni  is  a  unit  vector  normal  to  the  reaction  plane  and  along  the  direction 
defined  by  ki  x  k2.  If  now  such  a  beam  of  neutrons  is  elastically  scattered, 
the  differential  cross  section  of  the  scattered  neutrons  is  given  by: 

a(02,<}>)  =  a(02)[l  +  P1P2(ni  •  n2)] 


(1-1) 


t  I  "-  *  »  X-»«  «*  b8«u  Xino«o,[  3  ,M 
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target 


"  (f> 


Fig.  1-1.  Geometrical  arrangement  in  a  polarization  experiment 


where  a(02)  =  the  scattering  cross  section  if  the  incident  beam  of  neutrons 
were  unpolarized 

P2  =  polarization  of  the  scattered  neutrons  if  the  incident  beam 
of  neutrons  were  unpolarized 

n2  =  a  unit  vector  normal  to  the  scattering  plane  and  along  the 

f 

direction  k2  x  k2 
and  cos  <j>  =  nj  •  n2 

If  now  a  neutron  detector  is  placed  at  (02,O)  and  another  at  (02,tt),  the 
right-left  asymmetry,  e,  is  given  by 

a(02,O)  -  o(02,tt) 

£  ~  a  (02 ,0)  +  a  (02  ,tt) 

=  P1P2 


(1-2) 
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Since  e  can  be  measured  in  a  scattering  experiment  and  if  P2  is  known  also, 
the  polarization,  Pi  ,  of  the  neutrons  from  the  reaction  can  be  determined. 


1-2  A  Comparison  of  Different  Scatterers 


The  quantity,  P2 ,  is  frequently  referred  to  as  the  analysing  power  of 
the  scatterer.  One  of  the  methods  of  determining  P2  is  to  calculate  it 
from  phase  shifts  obtained  from  scattering  experiments  (Se  53,  Ho  66). 

Since  the  phase  shifts  are  best  known  for  ^He,  i2C  and  D0,  these  nuclei 
are  usually  used  as  scatterers.  Also  these  nuclei  have  zero  spin,  so  that 
for  a  given  incident  orbital  angular  momentum,  £,  there  are  only  two  partial 
wave  phase  shifts,  6  ,  corresponding  to  the  total  angular  momenta,  J  =  £  + 

Xj 

respectively.  The  analysing  power  can  then  be  written  as: 
p  =  21m (g*h) 

-  (1-3) 


where 


i£+ 

g  =  (1/k)  EPj^  (cos6)  [(£+l)sin  e  1 

£  1 


i<5£  , 

+  £  sin  6^  e  j 


(T)  +  -  1(6^  +  6£) 

h  =  (l/k)ZP^X0)sin(5^  -  6£)  e 


P  (cos0)  and  P^'(0)are  the  Legendre  and  associated  Legendre  polynomials  and  Im 

£  £ 

represents  the  imaginary  part  of  the  quantity  inside  the  brackets.  Hence  the 
analysing  power  of  the  scatterer  can  be  calculated  from  the  phase  shifts. 


Fig.  1-2  shows  the  analysing  powers  of  4He,  lZC  and  J^0  calculated  from 
the  known  phase  shifts  (Ha  63).  At  a  centre-of-mass  scattering  angle  of  135J, 


1 — 1 1  CM 
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Fig.  1-2. 


The  analysing  powers  of  4He,  12C  and  1 60  as  calculated 
from  the  known  phase  shifts  (Ha  63) 


En  (MeV) 
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the  analysing  power  for  uHe  is  large  and  slowly-varying  for  neutrons  of 
energy  from  3  to  20  MeV.  This  smooth  variation  with  energy  is  due  to  the 
two  very  broad  states,  p_  and  p  ,  in  5He  and  these  two  phase  shifts 

ill  !/2 

contribute  most  to  the  analysing  power  below  20  MeV.  The  analysing  power 
for  1  C  at  a  centre-of-mass  scattering  angle  of  45°  varies  drastically  at 
2.08  and  2.95  MeV.  This  is  due  to  the  two  narrow  states,  d5^2  and  d^2, 
in  i3C.  Between  2.2  and  2.8  MeV  the  analysing  power  is  large  and  slowly- 


varying,  so  l2C  is  a  good  analyser  in  this  energy  range.  As  for  iO0,  even 

for  neutrons  of  energy  below  1.5  MeV,  the  scattering  cross  section  shows 

three  well  separated  resonances  at  .43,  1.00  and  1.31  MeV  corresponding 
3—  2"*"  3— 

to  —  ,  —  and  —  states  in  1 70  respectively.  Consequently  the  analysing 
power  of  160,  as  shown  in  the  diagram  for  a  centre-of-mass  scattering  angle 
of  60°,  varies  quite  a  bit  over  this  energy  region.  Hence  of  the  three 
scatterers,  HHe  is  the  most  generally  useful  analyser. 

Another  advantage  of  using  4He  as  an  analyser  is  that  it  has  been 
found  to  scintillate  under  neutron  bombardment  (Th  59,  FI  59,  Si  61,  To  61, 
Es  65) .  Since  the  scintillation  signal  is  proportional  to  the  helium 
recoil  energy,  it  can  be  used  as  a  measure  of  the  neutron  bombarding  energy, 
Also  such  a  signal  can  be  used  to  set  up  a  time-of-f light  system,  which  is 
the  technique  adopted  in  our  system.  As  a  polarization  experiment  is 
essentially  like  a  double  scattering  experiment,  liquid  helium  is  used  in 


16, 


our  case  instead  of  helium  gas  in  order  to  improve  the  counting  rate.  Since 
liquid  helium  is  used,  a  specially  built  cryostat  must  be  used.  The  following 
section  describes  the  cryostat  employed  in  this  laboratory. 
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1-3  Liquid  Helium  Cryostat 

Fig.  1-3  shows  a  schematic  drawing  of  the  liquid  helium  cryostat, 
which  was  manufactured  by  Hofman  Cryogenics  (Si  61,  Ka  63).  It  is  a 
stainless  steel  dewar,  the  central  part  of  which  is  a  liquid  helium 
reservoir  of  1  litre  capacity.  Surrounding  it  is  a  liquid  air  reservoir 
to  reduce  the  temperature  gradient  between  the  liquid  helium  reservoir 
and  the  outside.  The  space  in  between  the  two  reservoirs  and  the  outer 
wall  is  under  high  vacuum.  A  copper  heat  shield  extends  from  the  bottom 
of  the  liquid  air  reservoir  down  to  the  thin  section  of  the  polarimeter. 

The  scattering  volume  of  the  polarimeter  is  at  the  bottom  end  of  the  liquid 
helium  reservoir,  with  a  wall  thickness  of  .010".  It  is  defined  by  a  .010" 
stainless  steel  cylinder  of  2"  in  diameter  and  4"  in  height.  A  hole  1/4" 
in  diameter  at  the  top  of  the  cylinder  allows  liquid  helium  to  pass  through. 

An  artificial  sapphire  window  2"  in  diameter  allows  the  scintillation 
signal  to  be  observed  by  a  RCA8575  phototube  placed  1/8"  beneath  the 
sapphire  window.  The  phototube  is  sealed  in  vacuum  by  means  of  an  0-ring 

1  —  7 

at  the  base  of  the  tube.  A  vacuum  of  about  2  x  10  '  mm  of  mercury  is  attained 
by  means  of  a  2"  Edwards  mercury  diffusion  pump  and  a  mechanical  rotary  fore- 
pump  . 

Since  the  scintillation  of  the  helium  recoil  is  in  the  ultraviolet  region 
(.1600  X),  a  wavelength  shifter  must  be  used  on  the  scattering  volume  to  bring 

Fig.  A2-1  in  Appendix  2  shows  a  schematic  drawing  of  the  mounting  of  the 
phototube. 
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Fig.  1-3 


A  schematic  drawing  of  the  liquid  helium  cryostat 
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the  signal  to  the  visible  region.  The  commonly  used  wavelength  shifters 

•j* 

are  pp ' -diphenylstilbene ,  quaterphenyl  and  popop.  pp ' -diphenyls tinlbene 
was  chosen  in  our  case,  since  it  was  found  to  give  the  best  energy  resolution 
when  tested  with  a  239Pu  alpha-particle  source  (Si  61).  In  order  to  improve 
the  light-collecting  properties  of  the  scattering  volume,  the  inside  of 
the  thin  stainless  steel  cylinder  was  first  coated  with  a  thin  layer  of 
aluminum  by  means  of  vacuum  evaporation,  followed  by  a  thin  protective  layer 
of  silicon  monoxide  and  then  a  fine  layer  of  white  paint.  The  paint  was 

ti 

made  by  mixing  one  part  of  zinc  oxide,  one  part  of  transparent  Tygon  paint 

ft 

and  two  parts  of  Tygon  thinner  .  Coarse  particles  were  removed  by  drain¬ 
ing  the  paint  through  a  funnel  filled  with  glass  wool.  The  paint  so  pre¬ 
pared  was  poured  onto  the  inside  of  the  cylinder  and  allowed  to  drain  and  dry 
before  the  next  layer  was  put  on.  Four  layers  of  paint  were  found  sufficient 
to  give  a  white  surface;  a  thick  layer  of  paint  was  found  to  crack  at  liquid 
helium  temperature.  The  wavelength  shifter  was  then  applied  by  slow  evapor¬ 
ation  to  a  thickness  of  about  100  yg/cm2.  The  sapphire  window  was  also 
lightly  coated  with  the  wavelength  shifter  to  a  thickness  of  about  50  yg/cm2. 
Since  the  pulse  height  of  the  helium  scintillations  was  found  to  deteriorate 
if  frozen  gases  or  water  vapour  were  collecting  at  the  bottom  of  the  scattering 
volume,  the  liquid  helium  reservoir  was  always  kept  under  an  atmosphere  of 
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helium  while  not  in  use.  Liquid  helium  was  transferred  at  the  top  of  the 
cryostat  by  means  of  a  double-walled  siphon.  A  brass  cone  was  attached 
on  top  of  the  1/4"  hole  of  the  thin  stainless  steel  cylinder  to  act  as 
a  siphon  guide  (see  Fig.  1-3)  so  that  one  arm  of  the  siphon  could  be 
inserted  into  the  scattering  volume  while  transferring,  as  this  would 
reduce  the  initial  splashing  of  liquid  helium. 

Two  carbon  resistors  of  16K  each  were  suspended  by  means  of  fine 
electric  wires  from  the  top  of  the  cryostat,  so  that  one  was  5"  and  the 
other  16"  from  the  bottom  of  the  scattering  volume.  These  two  resistors 
acted  as  a  liquid  helium  level  indicator,  as  their  resistance  increased 
to  2M  when  immersed  in  liquid  helium.  One  fill  of  liquid  helium  used 
about  3  litres  of  liquid  helium  and  lasted  for  about  3-1/2  hours.  The 
helium  gas  was  recovered  through  a  brass  pipe  leading  back  to  the  storage 
tank  of  the  liquefier. 

1-4  The  Analysing  Power  of  Helium 

Phase  shifts  for  n-^He  scattering  have  usually  been  derived  from  the 
more  accurately  known  p-4He  phase  shifts  by  accounting  for  the  difference 
in  the  Coulomb  energies  so  that  the  corresponding  energy  levels  in  5He  and 
5Li  are  matched.  The  Dodder-Gammel-Seagrave  phase  shifts  (Do  52,  Se  53)  for 
neutrons  of  energy  below  20  MeV  have  long  been  used  to  calculate  the  helium 
analysing  power.  Gammel,  Thaler  and  Perkins  determined  another  set  of 
phase  shifts  for  neutrons  of  energy  above  10  MeV  (Ga  58).  However,  these 
two  sets  of  phase  shifts  do  not  connect  smoothly  nor  do  they  fit  the  available 
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data  above  a  neutron  energy  of  15  MeV  (Ho  66) .  For  this  reason,  Hoop  and 
Barschall  (Ho  66)  determined  another  set  of  phase  shifts  covering  the 
neutron  energy  range  from  .5  to  30  MeV.  In  their  case,  the  n-4He  phase 
shifts  for  neutron  energy  above  17  MeV  were  derived  from  the  Weitkamp- 
Haeberli  p-4He  phase  shifts  (We  66)  and  the  s-and  p-wave  phase  shifts 
were  joined  smoothly  to  the  Dodder-Gammel-Seagrave  phase  shifts  at  low 
energies.  The  analysing  power  calculated  from  this  set  agrees  very  well 
with  polarization  measurements  using  the  polarized  neutrons  from  the 
T(p,n)3He,  D(d,n)3He  and  T(d,n)4He  reactions,  covering  the  neutron  energy 
range  6  to  24  MeV. 

For  neutrons  of  energy  below  8  MeV,  which  is  the  energy  range  of  our 
interest,  the  s-and  p-waves  are  most  important  in  the  partial  wave  analysis 
of  the  scattering  data,  as  the  d-waves  are  small  and  can  be  neglected.  Since 
the  s-and  p-wave  phase  shifts  can  be  determined  accurately,  the  calculated 
analysing  power  can  be  used  with  confidence  in  this  low  energy  range. 

At  a  neutron  energy  of  22.16  MeV,  the  17.669  MeV  state  of  5He  is  reached. 

3+ 

Since  J  =  —  for  this  state,  the  d phase  shift  varies  drastically  at  this 
energy.  Also  with  the  onset  of  the  4He(n,d)T  reaction,  the  phase  shifts 
become  complex  and  uncertain.  Hence  the  analysing  power  for  helium  was 
calculated  only  up  to  20  MeV  for  our  system.  The  Hoop-Bars chall  phase  shifts 
were  used,  and  since  the  highest  neutron  energy  used  was  20  MeV,  partial 
waves  above  d-wave  were  neglected.  For  convenience  in  the  computation  of  the 
analysing  power,  the  phase  shifts  were  expressed  as  a  power  series  in  neutron 
energy  by  means  of  a  least  squares  fitting  program^.  The  expressions  are 

^Written  by  W.  G.  Davies,  Ph.D.  thesis  (1966),  Department  of  Physics, 
University  of  Alberta. 
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given  below: 

1/2  3/2 

6 0  =  180  -  21.7684E  -  2.63613E  +  .68954E  ' 

for  E  _<  20  MeV 

8l  =  2.45763E  +  3.55942E2  -  .505754E3 

for  E  <_  2  MeV 

67  =  -10.9749  +  14.8534E  -  .97544E2  +  .0197815E3 

for  2  <  E  <_  20  MeV 

=  684.93  -  2312.49Ei/2  +  2493. 53E  -  807.582E3/2 

for  E  _<  2  MeV 

fit  =  .0648862  +  180.513E1/2  -  92.6748E  +  19.1394E3/2  -  1.43251E2 

for  2  <  E  <_  20  MeV 

67  =  -.0725766E  +  .0133187E2  for  8  <_  E  <_  20  MeV 

6|  =  . 128117E  -  . 0084296E2  +  .00110985E3 

for  8  <_  E  <_  20  MeV 

In  Fig.  1-4  the  points  represent  the  Hoop-Barschall  phase  shifts  while  the 
curves  represent  the  fitted  values  as  calculated  from  the  above  expressions. 

A  computer  program  has  been  written  to  calculate  the  n-^He  analysing 
power  for  a  finite  size  scatterer  and  a  finite  size  detector.  A  descrip¬ 
tion  and  a  listing  of  the  program  are  given  in  Appendix  1. 


n-4He  phase  shifts 
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CHAPTER  II 

EXPERIMENTAL  ARRANGEMENT 

2-1  Production  of  a  1  ns  Pulsed  Beam 

As  a  neutron  time-of-f light  technique  was  employed  in  the  polarization 

experiment,  a  brief  description  of  the  production  of  a  1  ns  pulsed  beam 

i 

with  the  5.5  MeV  Van  de  Graaff  accelerator  at  the  University  of  Alberta 
will  be  given. 

The  ionized  deuterium  gas  is  expelled  from  a  R.F.  ion  source  by 
applying  a  positive  voltage  to  the  probe  (see  Fig.  2-1).  After  focussing, 
the  ions  are  bent  through  an  angle  of  25 3  by  means  of  a  terminal  analysing 
magnet.  Two  sets  of  parallel  R.F.  deflection  plates  (1  MHZ)  perpendicular 
to  each  other  sweep  the  ion  beam  in  an  elliptical  orbit,  and  only  that 
portion  of  the  beam  sweeping  across  a  small  aperture  is  accepted.  In  this 
way  a  10  ns  pulsed  beam  is  produced.  The  pulsed  beam  is  then  accelerated 
down  the  column  of  the  accelerator.  A  90°  magnet  bends  the  beam  horizont¬ 
ally  and  deuterons  of  the  correct  energy  are  allowed  to  pass  through  a  pair 
of  energy  defining  slits.  After  passing  through  a  set  of  R.F.  deflection 

plates  (10  MHz) ,  a  beam  pulse  is  tilted  to  a  correct  angle  to  be  accepted 

ii 

by  a  Mobley  bunching  magnet  .  Since  the  front  part  of  the  pulse  traverses 

i 

Model  GN  Van  de  Graaff,  High  Voltage  Engineering  Corporation, 
Burlington,  Massachusetts, 
ft 
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Fig.  2-1.  Block  diagram  of  the  Van  de  Graaff  accelerator  and 
the  Mobley  magnet  system 
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a  longer  orbit  than  the  end  of  the  pulse,  the  10  ns  pulse  is  compressed 
into  a  1  ns  or  shorter  pulse  (Da  66A) .  The  1  ns  pulses  are  picked  off  by 
means  of  a  cylindrical  capacitor  placed  in  front  of  the  target  chamber 
and  are  used  as  stop  pulses  for  the  neutron  time-of-f light  system. 

The  average  beam  current  is  measured  at  the  insulated  target,  which 
is  connected  to  a  +300  volt  source,  to  prevent  secondary  electrons  pro¬ 
duced  at  the  target  from  escaping.  The  maximum  average  beam  current  used 
in  the  present  experiment  was  about  3  ya  corresponding  to  a  peak  current 
of  about  3  ma  for  a  repetition  rate  of  1  MHz.  Damage  to  the  target  would 

be  serious  if  a  higher  beam  current  were  used. 

2-2  Experimental  Configuration 

A  layout  of  the  experimental  set-up  is  shown  in  Fig.  2-2 a.  The  neutrons 

produced  at  the  target  were  first  scattered  by  the  scattering  volume  of 

the  polarimeter  and  then  detected  by  the  two  side  neutron  detectors.  The 

neutron  flight  path  between  the  target  and  the  polarimeter  was  set  at  2  m. 

The  two  neutron  detectors  were  each  set  at  an  angle  of  123°  (lab)  with 

respect  to  the  direction  of  the  incident  neutron  beam  and  the  distance  was 

25  cm  from  the  centre  of  the  polarimeter  to  the  front  face  of  the  neutron 

detector.  In  order  to  prevent  neutrons  produced  at  the  target  from  hitting 

the  neutron  detectors  directly,  two  lead  shadow  bars  each  60  cm  long  by 

20  cm  wide  by  25  cm  high  were  used.  The  polarimeter,  the  two  neutron 

detectors  and  the  two  shadow  bars  were  mounted  on  a  cart?  which  could  be 
i 

A  schematic  drawing  of  the  cart  is  shown  in  Fig.  A2-2  in  Appendix  2. 
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Fig.  2-2 a.  A  layout  of  the  experimental  arrangement 
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rotated  from  -20°  to  +120°  about  the  centre  of  the  target  chamber.  Also 
the  two  neutron  detectors  could  be  rotated  about  the  centre  of  the 
polarimeter  from  75°  to  130° „ 

The  two  neutron  detectors  were  identical  in  size  and  shape.  They 

1 

were  two  NE213  liquid  scintillators  of  the  bubble  free  type  ,  so  that 
they  could  be  placed  at  any  position  without  worrying  about  the  expansion 
bulb  on  the  side.  The  internal  dimensions  of  each  of  the  scintillators 
were  3-1/2"  in  diameter  and  2"  in  height.  They  were  mounted  on  two 
Philips  5"  XP1040  phototubes.  As  well  as  the  y-metal  shield  around  the 
phototube  provided  by  the  manufacturer,  another  thick  y-metal  shield  was 
put  around  the  assembly  of  the  scintillator  and  the  phototube  to  reduce 
to  a  minimum  the  gain  shift  of  the  phototube  due  to  stray  magnetic  fields. 
This  thick  y-metal  shield  served  as  a  light  shield  as  well. 

Two  pictures  of  the  polarization  system  taken  at  two  different  angles 
are  shown  in  Fig.  2-2b.  The  experimental  configuration  can  be  seen  quite 
clearly  from  these  pictures. 

2-3  Electronic  Configuration 

Fig.  2-3  shows  a  block  diagram  of  the  electronics  of  the  polarization 
system.  A  neutron  time-of-f light  technique  is  used  to  resolve  neutrons 
of  different  energies.  The  scintillation  signal  from  the  scattering  volume 
of  the  polarimeter  is  used  as  a  start  signal  for  the  time-of-f light  system, 
while  the  pulse  from  the  Mobley  capacitor  pick-off  is  used  as  a  stop  signal. 
The  electronics  for  the  two  side  neutron  detectors  are  identical.  A  fast 
coincidence  is  set  between  the  polarimeter  and  each  of  the  two  side  detectors 

^Nuclear  Enterprises  Corporation,  Winnipeg,  Manitoba 
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Fig,  2-2b,  Pictures  of  the  polarization  system 

The  top  picture  shows  the  whole  assembly  as  viewed  from  the 
direction  of  the  target  chamber.  At  the  bottom  of  the 
picture  is  the  stainless  steel  target  chamber;  beyond  that 
are  the  two  neutron  detectors  marked  A  and  B.  The  two  lead 
shadow  bars  are  in  between  the  target  chamber  and  the  neutron 
detectors.  The  liquid  helium  polarimeter  is  seen  straight 
ahead  of  the  target  chamber.  The  lower  picture  shows  a  side 
view  of  the  assembly.  In  this  picture  the  mercury  diffusion 
pump  and  the  rotary  fore-pump  are  seen  on  the  left-hand  side. 
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Fig.  2- 


A  block  diagram  of  the  electronic  arrangement  of 
the  polarization  system 
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for  the  measurement  of  right-left  asymmetry.  The  output  of  dynode  13  from 
each  of  the  detectors  is  used  to  set  the  neutron  threshold.  A  detailed 
description  of  the  system  is  given  in  the  following  paragraphs. 

The  time-to-amplitude  conversion  system  was  designed  by  D.  L.  Wieber 
(Wi63).  Fig.  2-4  shows  a  block  diagram  of  the  system.  A  tunnel  diode 


Start  Pulse 


Fig.  2-4.  Block  diagram  of  the  time-to-amplitude  conversion  system. 

trigger  is  used  at  the  input  stage  of  both  the  start  and  stop  discriminators, 
whose  thresholds  can  be  set  from  10  to  600  mV.  They  are  stable  up  to  a 
repetition  rate  of  10  MHz.  The  output  pulse  from  the  start  discriminator 
switches  on  a  constant  current,  I,  to  charge  a  mica  capacitor,  C,  which  is 
kept  charged  until  a  pulse  from  the  stop  discriminator  comes  along  to  switch 
it  off.  If  the  time  interval  between  the  start  and  stop  pulses  is  t,  the 
voltage,  Vc,  across  the  capacitor  will  be  Vc  =  t.  Since  I/C  remains 
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constant,  a  time  interval  can  be  converted  linearly  into  a  related  pulse 
height.  The  integral  linearity  of  the  system  was  found  to  be  about  1%. 

In  the  polarization  experiment,  a  RCA8575  phototube  is  used  to  observe 
the  scintillation  signals  in  the  scattering  volume  of  the  polarimeter  and 
the  anode  pulses  are  used  as  start  signals  for  the  time-of-f light  system. 
Even  for  a  monoenergetic  source  of  neutrons  bombarding  the  polarimeter,  the 
pulse  height  would  still  vary  with  the  direction  of  the  helium  recoil. 

Hence  a  time  spread  would  result  as  different  size  pulses  cross  the  thres¬ 
hold  of  the  start  discriminator,  even  though  the  threshold  could  be  set  as 
low  as  30  mV.  This  anode  time  shift  is  the  chief  limitation  for  the  time 
resolution  of  the  system. 

Dr.  W.  K.  Dawson  of  this  laboratory  has  written  a  general  purpose 
kicksorter  program  (GPKS)  for  the  on-line  use  of  the  SDS  920  computer. 

There  are  two  routines,  called  TWOD  and  TOFA,  in  it  for  the  compensation 
of  this  anode  time  shift.  The  TWOD  routine  is  to  generate  a  shift  table 
to  be  used  for  the  TOFA  routine,  which  gives  a  corrected  time-of-f light 
spectrum  for  neutrons  of  different  energy  groups.  A  description  of  these 
two  routines  will  be  given  in  the  next  paragraph. 

Fig.  2-5  shows  the  underlying  principle  in  the  TWOD  routine.  Pulses 
a,  b  and  c  represent  events  due  to  neutrons  of  the  same  energy  but  with 
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different  helium  recoil  directions.  They  should  all  be  detected  at  the 
same  time.  Of  these  three  pulses,  pulse  a  crosses  the  discriminator 
threshold  at  an  earlier  time  than  both  b  and  c.  Consequently  on  the  time- 
of-flight  spectrum,  events  b  and  c  would  be  treated  as  neutrons  of  lower 
energy  and  the  peak  corresponding  to  this  neutron  group  would  be  quite 
broad.  In  order  to  compensate  for  this  anode  time  shift,  the  corresponding 
pulses  from  a  dynode  output  are  examined.  A  two  parameter  correlation 
between  the  anode  (time)  pulse  and  the  dynode  (amplitude)  pulse  is  recorded 
to  determine  the  time  shift  as  a  function  of  pulse  amplitude.  In  this  way 
a  shift  table  is  generated  to  be  used  in  the  TOFA  routine.  In  this  routine 
the  flight  time  of  each  event  is  shifted  the  right  amount  to  give  a  corrected 
time-of-f light  spectrum  according  to  the  amplitude  of  the  dynode  pulses. 

A  monoenergetic  source  of  neutrons  from  the  reaction  T(d,n)4He  was  used 
to  generate  a  shift  table  for  the  polarimeter.  The  deuteron  bombarding 
energy  was  3  MeV  and  the  polarimeter  was  set  at  0°  with  respect  to  the 
direction  of  the  incident  deuterons,  so  that  the  neutron  energy  was  19.58 
MeV.  The  time-of-f light  system  for  the  neutrons  was  set  between  the  target 
and  the  polarimeter  with  a  flight  path  of  2m.  A  block  diagram  of  the 
electronics  employed  is  shown  in  Fig.  2-6.  The  anode  output  from  the 
RCA8575  photomultiplier  tube  was  used  as  the  start  pulse  while  the  Mobley 
magnet  capacitor  pick-off  was  used  as  the  stop  pulse  for  the  time-of-f light 
system.  The  time-to-amplitude  converter  output  was  fed  into  ADC(A)  after 
going  through  a  delay  line  amplifier.  The  11*-^  dynode  output  from  the 
photomultiplier  was  fed  into  a  double  delay  line  amplifier  and  then  into 
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Fig.  2-6.  An  electronic  block  diagram  for  the  generation  of  a 
shift  table 
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ADC(B).  The  double  delay  line  amplifier  output  was  also  fed  into  a  discrim¬ 
inator  so  that  a  lower  threshold  for  the  polarimeter  could  be  set  if  desired. 
The  discriminator  output  was  then  fed  into  a  coincidence  interrupt,  which 
in  turn  was  used  to  gate  the  two  ADC’s.  The  purpose  of  using  the  coincidence 
interrupt  was  to  make  sure  that  both  ADC's  and  the  computer  would  all 
accept  or  reject  a  particular  event. 

The  range  for  ADC(B)  was  1024  channels,  and  these  1024  channels  were 
divided  into  56  (amplitude)  bins  for  the  shift  correction.  Since  a  smaller 
pulse  crosses  the  start  discriminator  threshold  with  a  smaller  slope  than  a 
larger  one,  the  time  spread  is  larger  among  the  smaller  amplitude  group. 

For  this  reason  the  bin  size  for  the  lower  bins  was  smaller  than  those  for 
the  upper  bins,  so  that  the  time  resolution  of  the  neutron  peak  in  the  lower 
bins  would  not  be  too  poor.  Starting  from  bin  0,  the  first  32  bins  were  of 
8  channels  wide,  the  next  16  bins  of  16  channels  wide,  and  the  next  8  bins 
of  64  channels  wide.  In  this  way  the  time-of-f light  spectrum  for  this  mono- 
energetic  neutron  group  was  divided  into  56  bins.  As  can  be  seen  from 
Fig.  2-7  which  shows  a  picture  of  the  TWOD  spectrum  taken  directly  from 
the  display  unit  of  the  on-line  computer  system,  the  peak  positions  of  the 
different  bins  shifted  to  lower  channels  with  respect  to  the  highest  bin. 
Hence  in  constructing  a  shift  table,  the  peak  position  of  the  highest  bin 
was  taken  as  the  reference  and  the  other  peak  positions  were  recorded  with 
respect  to  this  reference  position.  In  order  that  all  the  peak  positions 
would  occur  at  the  same  channel  as  the  highest  bin,  the  time-of-f light 
spectrum  in  different  bins  was  incremented  by  a  number  of  channels  equal  to 
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Fig.  2-7.  A  picture  of  the  TWOD  spectrum  taken  directly  from  the  display 
unit  of  the  on-line  computer  system. 

The  time  axis  is  along  the  X-axis  with  decrease  in  flight  time 
along  the  positive  direction  of  the  X-axis.  The  amplitude  axis 
is  along  the  Y-axis. 
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the  difference  in  peak  positions  between  the  reference  peak  and  the  peak 
in  the  bin.  The  shift  table  generated  for  the  polarimeter  is  shown  in 
Table  2-1  and  a  plot  of  this  shift  table  is  shown  in  Fig.  2-8.  From  the 
TWOD  spectrum,  the  time  resolution  for  the  different  bins  was  extracted. 

In  Fig.  2-9,  the  time  resolution  has  been  converted  from  channels  to  nano¬ 
seconds  and  the  bins  have  been  converted  into  helium  recoil  energy  (also 
see  Table  2-1) .  The  time  resolution  for  helium  recoil  energy  below  1  MeV 
was  quite  poor,  while  for  the  recoil  energy  range  from  1  to  12  MeV,  the 
time  resolution  improved  gradually  from  about  2  ns  to  1  ns.  Fig.  2-10 
shows  a  time-of-f light  spectrum  of  the  19.58  MeV  neutron  group  from  the 
T(d,n)l+He  reaction.  Without  the  anode  time  shift  correction,  the  time 
resolution  (FWHM)  was  4.7  ns,  while  that  with  a  time  shift  correction 
was  lo61  ns.  The  polarimeter  threshold  was  set  to  accept  pulses  due  to 
helium  recoils  of  .5  MeV  or  higher  in  energy;  this  corresponds  to  an 
incident  neutron  energy  of  about  1  MeV.  For  best  time  resolution,  it  was 
found  that  the  start  discriminator  had  to  be  set  just  above  the  noise  level 
of  the  phototube,  unlike  the  case  of  a  fast  liquid  scintillator,  like  NE213. 

Since  liquid  helium  also  scintillates  with  y-rays,  the  . 511y  from 
22Na  was  used  to  investigate  the  time  resolution  of  the  polarimeter  for 
helium  recoil  energy  below  ,5  MeV.  A  22Na  y-ray  source  was  put  between  the 
polarimeter  and  a  small  NATON  organic  phosphor  mounted  on  a  Philips  56AVP 
phototube.  The  electronic  arrangement  was  similar  to  that  shown  in  Fig.  2-6 
except  that  the  anode  signal  from  the  56  AVP  was  the  source  of  the  stop 
pulse  instead  of  the  Mobley  capacitor  pick-off  and  a  narrow  window  was  set 
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Table  2-1. 

Shift  Table 

and  Time  Resolution 

for  the  Polarimeter 

Bin  No. 

Time  Shift 
(Channels) 

a-Recoil  Energy 
(MeV) 

Time  Resolution 
(ns) 

0 

62 

.049 

too  wide 

1 

57 

.147 

too  wide 

2 

52 

.245 

4.36 

3 

47 

.343 

4.18 

4 

42 

.441 

3.90 

5 

39 

.538 

3.45 

6 

35 

.  636 

2.91 

7 

33 

.734 

2.63 

8 

31 

.832 

2.27 

9 

30 

.930 

2.09 

10 

29 

1.028 

2.09 

11 

27 

1.126 

1.91 

12 

27 

1.224 

1.82 

13 

25 

1.322 

1.82 

14 

25 

1.420 

1.54 

15 

24 

1.517 

1.63 

16 

23 

1.615 

1.63 

17 

23 

1. 713 

1.54 

18 

23 

1.811 

1.45 

19 

22 

1.909 

1.45 

20 

21 

2.007 

1.45 

21 

21 

2.105 

1.54 

22 

20 

2.203 

1.36 

23 

20 

2.301 

1.54  i 

24 

19 

2.399 

1.27  1 

25 

19 

2.496 

1.45 

26 

19 

2.594 

1.18 

.  - 

Table  2-1  (continued) 
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Bin  No. 

Time  Shift 
(Channels ) 

a-Recoil  Energy 
(MeV) 

Time  Resolution 
(ns) 

27 

18 

2.692 

1.27 

28 

18 

2.790 

1.45 

29 

17 

2.888 

1.36 

30 

17 

2.986 

1.36 

31 

17 

3.083 

1.27 

32 

17 

3.230 

1.18 

33 

16 

3.426 

1.27 

34 

16 

3.622 

1.27 

35 

15 

3.818 

1.27 

36 

14 

4.014 

1.27 

37 

14 

4.209 

1.27 

38 

14 

4.405 

1.27 

39 

13 

4.601 

1.18 

40 

13 

4.797 

1.18 

41 

12 

4.993 

1.18 

42 

12 

5.188 

1.09 

43 

12 

5.384 

1.27 

44 

12 

5.580 

1.18 

45 

11 

5.776 

1.27 

46 

11 

5.972 

1,18 

47 

11 

6.167 

1.09 

48 

10 

6. 656 

1.18 

49 

9 

7.439 

1.09 

50 

8 

8.223 

1.09 

51 

8 

9.006 

1.09 

52 

7 

9.789 

1.00 

53 

5 

10.572 

1.09 

54 

5 

11.355 

1.09 

55 

0 

12.138 

1.00 

Fig.  2-8. 
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Time  resolution  of  the  polarimeter  against  helium  recoil 
energies,  using  the  19 <,58  MeV  neutrons  obtained  from  the 
T(d,n)4He  reaction  at  Ej  =  3.0  MeV  and  0n  =  0° 


Fig.  2-9. 
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Corrected  and  uncorrected  neutron  time-of-flight  spectra  from 
the  T(d,n)4He  reaction  with  =  3  MeV,  0n  =  0°  and  a  flight 
path  of  2  metres. 


Fig.  2-10. 
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on  the  13^  dynode  output  of  the  56  AVP  to  gate  the  ADC's,  so  that  the  time 
resolution  of  the  system  would  be  mostly  due  to  the  polarimeter.  The  time 
resolution  is  shown  in  Fig.  2-11.  In  order  to  express  the  abscissa  in 
terms  of  helium  recoil  energy,  the  pulse  height  of  the  .511  MeV  y-ray  and 
that  due  to  5.1  MeV  229Pu  a-particles  were  compared. 

Although  the  energy  resolution  of  the  system  was  achieved  entirely  by 
the  time-of-f light  technique,  the  pulse  height  resolution  from  a  linear 
channel  was  also  investigated.  A  239Pu  a-source  was  introduced  into  the 
scattering  volume  of  the  polarimeter  by  means  of  a  long  stainless  steel 
rod  with  a  teflon  section  to  reduce  heat  conduction.  The  energy  resolution 
(AE/E)  was  found  to  range  between  25  and  35  %  depending  on  the  position  of 
the  a-source.  Fig.  2-12  shows  a  plot  of  the  energy  resolution  against  the 
a-source  position.  The  energy  resolution  deteriorated  as  the  distance 
between  a-source  and  window  was  increased.  This  indicates  that  the  light 
collecting  properties  of  the  stainless  steel  cylinder  are  poor.  Visually, 
the  sapphire  window  seemed  to  be  coated  rather  too  thickly  with  the  wave¬ 
length  shifter.  This  affects  the  energy  resolution,  since  considerable 
light  from  the  scintillation  event  is  reflected  back  into  the  scattering 
volume.  It  is  possible  that  improving  these  two  things  might  improve  the 
time  resolution  as  well.  Since  the  only  access  to  scattering  volume  is  by 
unsoldering  the  tail  section  of  the  polarimeter,  and  it  has  been  found  to 
be  difficult  to  resolder  it  leak  tight,  the  reflector  has  not  been  changed 


since  it  was  installed. 
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Fig,  2-11,  Time  resolution  of  the  polarimeter  using  the  ,511  MeV 
y-rays  of  22Na 
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Fig.  2-12.  Energy  resolution  of  the  polarimeter  obtained  from  a  linear 
channel  against  the  position  of  a  239Pu  a-particle  source 
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A  detailed  description  of  the  electronic  block  diagram  of  the  polari¬ 
zation  system  will  now  be  given  (see  Fig.  2-3).  The  anode  output  from  the 
RCA8575  phototube  was  fed  into  the  start  discriminator  while  the  Mobley 
magnet  capacitor  pick-off  was  first  fed  into  a  fast  amplifier  system  (Hewlett 
Packard  461A  and  462A)  and  then  to  the  stop  discriminator  of  the  time-of- 
flight  system.  The  time  range  of  the  time-to-amplitude  converter  was  300 
nano-seconds.  A  delay  line  amplifier  was  used  to  convert  the  negative 
time-to-amplitude  converter  output  to  a  bipolar  pulse  so  that  the  pulse 
shape  would  be  of  the  required  shape  to  be  accepted  by  ADC (A).  The  11^ 
dynode  output  from  the  RCA8575  phototube  was  first  fed  into  a  double  delay 
line  amplifier  and  then  into  ADC(B).  The  outputs  of  both  ADC(A)  and  ADC(B) 
were  fed  into  the  computer  for  the  on-line  correction  of  the  anode  time 
shift.  All  these  formed  the  main  time-of-f light  system  to  resolve  neutrons 
of  different  energy  groups. 

Besides  feeding  into  ADC(B),  the  output  of  the  double  delay  line  ampli¬ 
fier  was  also  fed  into  a  linear  amplifier  and  then  to  a  discriminator,  so 
that  a  neutron  threshold  could  be  set  for  the  polarimeter.  The  purpose  of 
using  the  linear  amplifier  was  to  achieve  a  lower  threshold  than  that 
obtainable  by  feeding  the  double  delay  line  amplifier  output  directly  to 
the  discriminator.  The  threshold  set  for  the  polarimeter  corresponded  to 
about  250  keV  of  helium  recoil  energy. 

The  electronic  systems  for  the  two  neutron  detectors  were  identical. 

The  anode  output  of  the  XP1040  phototube  was  doubly  RC  differentiated,  so 
that  the  amplitudes  of  the  positive  and  the  negative  portions  of  the  pulse 
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were  about  equal  and  of  about  10  nanosecond  width.  The  pulse  was  then 
inverted  by  a  fast  pulse  transformer  ,  so  that  the  cross  over  of  the 
pulse  was  used  instead  of  the  leading  edge  to  trigger  the  start  discrim¬ 
inator,  since  there  was  no  anode  time  shift  compensation  done  for  the 
neutron  detectors  and  using  the  cross  over  would  result  in  a  smaller  anode 
time  shift  than  using  the  leading  edge.  The  stop  signal  for  the  fast 
coincidence  was  taken  from  the  fast  fan  out  of  the  polarimeter.  The 
time-to-amplitude  converter  output  was  first  fed  into  a  delay  line  ampli¬ 
fier  to  obtain  the  correct  pulse  shape  as  required  by  the  dual  channel 
analyser.  A  window  50  nanoseconds  wide  was  set  to  include  all  the 
neutron  peaks  and  the  y-ray  peak  by  means  of  the  dual  channel  analyser, 
which  also  provided  an  adjacent  upper  window  50  nanoseconds  wide  to  assess 
the  time  uncorrelated  background.  The  13^  dynode  output  from  the  XP1040 
was  first  fed  into  a  double  delay  line  amplifier  and  then  to  a  discrimin¬ 
ator  so  that  the  threshold  of  the  two  neutron  detectors  could  be  matched. 

Four  3-fold  coincidence  units  were  used,  so  that  the  true  and  the 
random  time-of-f light  spectra  for  the  neutrons  scattered  to  the  left  and 
the  right  could  be  recorded  all  at  one  time.  To  3-fold  coincidence  unit  A, 

the  discriminator  output  for  setting  the  threshold  of  the  polarimeter,  the 

i*f* 

true  fast  coincidence  for  neutrons  scattered  to  the  right  and  the 
f 

E  H  Research  Laboratory,  Inc.  ZT9-9X  pulse  transformer 
ft 

Since  the  Basel  sign  convention  is  used,  the  detector  on  the 
left  hand  side  of  the  page  is  actually  the  right  detector  according  to 
the  sign  convention. 
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discriminator  output  for  setting  a  threshold  for  the  right  neutron  detector 
were  fed.  For  3-fold  coincidence  unit  B,  two  of  the  inputs  were  the  same 
as  in  unit  A  and  the  third  one  was  from  the  random  fast  coincidence  windows 
for  the  right  detector.  3-fold  coincidence  units  C  and  D  played  the  same 
roles  for  the  left  neutron  detector.  The  outputs  of  the  four  3-fold 
coincidence  units  were  used  to  route  ADC(A),  so  that  the  first  and  the 
second  1024  channels  would  correspond  to  the  true  and  random  time-of-f light 
spectra  of  neutrons  scattered  to  the  right,  while  the  third  and  the  fourth 
1024  channels  would  correspond  to  the  respective  spectra  of  neutrons 
scattered  to  the  left.  Fig.  2-13  shows  a  picture  of  the  four  time-of-f light 
spectra  taken  directly  from  the  computer  display  unit  for  the  reaction 
28Si(d,n)29P. 

2-4  False  Asymmetry 

As  the  polarization  of  the  neutrons  from  a  reaction  is  measured  in  terms 
of  the  counts  recorded  by  the  right  and  the  left  detectors,  any  false 
asymmetry  due  to  defects  in  the  system  would  result  in  wrong  values  for 
the  polarization.  The  two  main  sources  of  false  asymmetry  are  called  geo¬ 
metric  asymmetry  and  electronic  asymmetry.  A  thorough  discussion  on  this 
subject  was  given  by  Haeberli  in  Chapter  V.G.  of  Fast  Neutron  Physics,  Part  II, 
and  the  discussion  below  will  be  entirely  concerned  with  the  elimination  of 


false  asymmetry  in  our  system. 
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Fig.  2-13.  A  picture  of  the  z8Si(d,n)29P  neutron  time-of-f light  spectra 
taken  directly  from  the  display  unit  of  the  on-line  computer 
system. 

The  picture  shows  the  four  1024  channels  corresponding  to  the 
true  and  random  time-of-f light  spectra  of  neutrons  scattered 
to  the  right  and  the  left  by  the  polarimeter. 

See  Fig.  5-2  on  page  70  for  the  identification  of  the  differ¬ 
ent  peaks o 
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(a)  Geometrical  Alignment 

The  whole  assembly  of  the  polarimeter,  the  two  neutron  detectors  and 
the  two  lead  shadow  bars  was  mounted  on  a  channel  iron  steel  cart  (see 
Fig.  2-2  and  Fig.  A2-2) .  The  centre  about  which  the  two  neutron  detectors 
were  rotated  was  adjusted  so  that  it  was  vertically  below  the  centre  of 
the  scattering  volume  of  the  polarimeter.  This  was  achieved  by  viewing 
the  centre  of  the  scattering  volume  by  two  transits  set  at  about  90°  apart. 
Adjustments  were  made  on  the  position  of  the  polarimeter  until  the  two 
centres  were  in  the  same  vertical  line  when  the  transits  were  tilted  up  and 
down  to  view  the  two  centres „ 

In  order  to  find  the  zero  degree  position  of  the  polarimeter  with 
respect  to  the  incident  beam  of  deuterons,  an  angular  distribution  of  the 
neutrons  from  the  DCdjnj^He  reaction  was  measured  from  -20°  to  +20°  at  5 
MeV  bombarding  energy.  Since  the  neutron  angular  distribution  peaks  very 
strongly  at  0°,  the  0°  position  for  the  polarimeter  was  determined  accurately 
from  the  peak  of  the  angular  distribution.  The  0°  position  was  found  to 
reproduce  with  +  1°. 

(b)  Electronic  Alignment 

The  .511  and  the  1.28  MeV  y-rays  from  2^Na  were  used  to  set  the  thres¬ 
holds  and  gains  of  the  two  neutron  detectors. 

Fig.  2-14  shows  a  block  diagram  of  the  electronics.  The  shape  of  the 
y-ray  spectrum  is  shown  in  Fig.  2-15. 
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22Na 
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Fig.  2-14.  A  block  diagram  for  setting  the  threshold  of  the  neutron 
detectors . 
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Fig.  2-15.  z2Na  y-ray  spectrum 


SC  ,.=•  I  .  ►  'M  -l  3i  1 


. 


■ 


■ 

43 

The  discriminator  level  was  adjusted  to  obtain  the  value  of  Cco  so  that 
the  following  ratio  was  obtained 

C  -  C 

D  _  1-28  .511  o  o 

K  _  p  3  •  Z 

*511  CO 

With  such  a  ratio,  the  threshold  would  correspond  to  about  .5  MeV  of  neutron 
energy . 

The  .511  MeV  y-rays  from  ^2Na  were  also  used  in  setting  up  the  50  nano¬ 
second  fast  coincidence  windows.  A  time-of-f light  spectrum  was  recorded 
by  putting  the  y-ray  source  between  the  polarimeter  and  one  of  the  neutron 
detectors.  The  peak  corresponding  to  the  .511  MeV  y-rays  gives  essentially 
zero  time,  since  y-rays  travel  much  faster  than  any  neutron  group  of  interest. 
Since  the  distance  between  the  polarimeter  and  the  neutron  detector  was  only 
25  cm,  the  time-of-f light  spectrum  for  neutrons  of  a  few  MeV  overlapped  with 
y-ray  peak.  For  this  reason  the  lower  window  of  the  dual  channel  analyser 
for  the  true  events  was  set  50  nanoseconds  wide  to  include  the  upper  tail 
of  the  y-ray  peak,  while  the  lower  edge  of  the  upper  window  was  set  at  the 
upper  tail  of  the  y-ray  peak  to  cover  a  region  of  50  nanoseconds  for  the 
assessment  of  the  time  uncorrelated  random  background. 

Even  after  the  above  adjustments,  it  was  usually  found  that  there  was 
some  false  asymmetry.  In  order  to  eliminate  this,  the  two  neutron  detectors 
and  the  electronics  associated  with  each  of  them  were  interchanged  in  their 
left  and  right  positions.  The  underlying  principle  is  illustrated  in  the 


following  paragraph. 
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Suppose  the  two  neutron  detectors  are  marked  A  and  B.  Let  R  and  L 
be  the  numbers  of  neutrons  scattered  to  the  right  and  the  left  by  the 
polarimeter  and  detected  by  the  right  and  the  left  detectors  respectively, 
then  ideally  the  right-left  ratio,  r,  is 


r  =  R/L  . 


If,  however,  the  efficiencies,  which  are  affected  by  both  the  detector 
threshold  and  the  fast  coincidence  window,  of  the  detectors  are  not 
exactly  the  same,  R/L  is  not  equal  to  r.  Let  the  number  of  neutrons 
detected  by  detector  A  when  situated  at  right  be  R^  and  that  by  detector 
B  situated  at  left  be  Lg,  then  instead  of  the  true  right-left  ratio,  one 
measures 


ra/Lb  =  r6 


where  6  is  due  to  the  difference  in  efficiencies.  If  now  the  two  detectors 
are  interchanged  in  their  positions,  the  counts  recorded  by  them  will  be 
Rg  and  L^  and  we  have 


rb/la  =  r/6 


Hence  the  true  right-left  ratio  is  obtained  if  the  geometrical  mean  of  the 


observed  right-left  ratios  is  taken,  i.e. 


r 


/ 


ra  rb 
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lb  la 


This  is  a  very  powerful  way  of  eliminating  not  only  the  electronic  asymmetry, 
but  also  the  false  asymmetry  due  to  the  two  detectors  not  being  placed  with 
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perfect  symmetry  about  the  polarimeter,  provided  the  right  and  the  left 
positions  of  the  detectors  remain  fixed. 

Instead  of  changing  the  two  detectors,  the  same  effect  is  attained 
if  the  polarimeter  angle  is  changed  from  0  to  -9,  since  the  two  detectors 
change  role  according  to  the  sign  convention.  For  0  equal  or  less  than  20°, 
measurements  were  done  in  this  way  instead  of  moving  the  detectors.  Measure¬ 
ments  from  both  methods  agreed  to  within  2  %  in  most  cases. 

2-5  The  Effectiveness  of  the  Shadow  Bars 

In  order  to  investigate  the  effectiveness  of  the  lead  shadow  bars 
(60  cm  long  by  20  cm  wide  by  25  cm  high) ,  neutrons  from  the  reaction 
T(d,n)4He  at  a  deuteron  bombarding  energy  of  3.5  MeV  were  used.  The  polar¬ 
imeter  was  set  at  0°  with  respect  to  the  direction  of  the  incident  beam  of 
deuterons  and  at  a  distance  of  1.5  metres  away  from  the  target.  A  time-of- 
flight  spectrum  was  recorded  between  the  polarimeter  and  one  of  the  neutron 
detectors,  which  was  set  at  90°  with  respect  to  the  direction  of  the  incident 
neutrons  and  at  a  distance  of  55  cm  away  from  the  polarimeter.  Time-of-f light 
spectra  for  neutrons  with  and  without  the  shadow  bar  are  shown  in  Fig.  2-15. 
Without  the  shadow  bar  the  background  to  signal  ratio  was  found  to  be  25.5  %, 
while  with  the  shadow  bar  between  the  target  and  the  neutron  detector,  this 
ratio  was  reduced  to  3.9  %.  It  was  also  noticed  that  in  the  presence  of  the 
shadow  bar,  the  width  of  the  peak  (FWHM)  improved  from  6.3  to  5.8  nanoseconds. 
It  should  be  mentioned  that  there  was  no  anode  time  shift  compensation  for  the 
time-of-f light  system  and  the  time  uncertainty  at  which  the  neutrons  were 
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Fig.  2-16.  Time-flight-spectra  to  show  the  effectiveness  of  the  shadow 
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detected  was  roughly  increased  by  a  factor  of  two  over  that  of  a  conven¬ 
tional  time-of-f light  experiment,  as  the  neutrons  had  to  traverse  the 
scattering  volume  of  the  polarimeter  to  give  a  start  signal  as  well  as  the 
neutron  detector  to  give  a  stop  signal  for  the  time-of-f light  system. 
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CHAPTER  III 

CALIBRATION  OF  THE  POLARIZATION  SYSTEM  WITH 
D(d,n)3He  NEUTRONS 

The  study  of  the  polarization  of  neutrons  from  the  D(d,n)3He  reaction 
has  long  been  of  interest.  Blin-Stoyle  (B1  51)  analysed  the  reaction 
by  the  method  of  partial  waves  and  arrived  at  the  following  expression 
for  the  prediction  of  polarization  for  the  outgoing  particles. in  a  nuclear 
reaction: 

L 

max  q) 

P(0)a(9)  =  l  Al  Pl  (cos6)  (3-1) 

L=1 

where  P(6)  =  polarization  of  the  outgoing  particles  emitted  at  an  angle  0 

a(0)  =  differential  cross  section 

Lmax  =  highest  partial  wave  used  in  the  reaction 

A^  =  a  coefficient  depending  on  the  matrix  elements  describing 
the  reaction 

(1) 

and  Pl  (cos0)  =  associated  Legendre  polynomial. 

If  the  largest  effective  values  of  the  incident  orbital  angular  momentum, 
final  orbit  angular  momentum  and  total  angular  momentum  are  £,  £’  and  J, 
then  Lmax  in  equation  (3-1)  must  be  less  than  or  equal  to  2£,  2£*  and  2J. 
In  addition,  if  the  interfering  states  have  the  same  parity,  all  the 


L- values  must  be  even. 
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In  the  D(d,n)3He  reaction,  since  both  the  incident  particles  and  the 


target  nuclei  are  identical  and  only  even  L-values  contribute,  equation 
(3-1)  can  be  written  as: 


(3-2) 


n 


The  determination  of  the  coefficients  an  gives  the  strength  of  the  spin- 
orbit  coupling  in  the  reaction.  Apart  from  the  determination  of  these 
coefficients,  polarization  experiments  have  been  performed  to  investigate 
the  possibility  of  using  this  mono-energetic  neutron  source  as  a  polarized 
neutron  source. 

In  1961  Dubbeldam  and  Walter  (Du  61)  made  a  survey  of  the  polariza¬ 
tion  measurements  from  this  reaction.  As  seen  from  equation  (3-2) , 
P(0)a(0)  is  largest  for  0  =  45°(c.m.).  For  this  reason  most  of  the 
polarization  measurements  were  performed  near  this  angle.  Dubbeldam  and 
Walter  found  that  for  deuteron  energies  above  2  MeV  their  measurements 
at  45°(c.m.)  were  in  disagreement  with  the  previous  measurements  around 
this  angle.  Since  that  time,  more  data  have  been  obtained  in  different 
laboratories.  In  1965  Purser,  Sawers  and  Walter  (Pu  65)  made  another 
survey  of  the  situation.  A  summary  of  their  survey  is  given  in  the  next 
paragraph. 

Fig.  3-1  presents  the  results  of  most  of  the  polarization  measure¬ 
ments  of  neutrons  from  the  D(d,n)3He  reaction  around  45°(c.m.).  Below 
deuteron  energies  of  2  MeV,  all  measurements  agree  quite  well  (Me  54, 

Le  57,  Pa  58),  while  above  2  MeV,  the  measurements  separate  into  two 
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Fig,  3-1  Polarization  measurements  of  neutrons  from  the  D(d,n)  He 
reaction  around  450(c.m„)° 
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groups.  The  first  reported  polarization  measurement  for  this  reaction 
above  2  MeV  was  performed  by  Daehnick  (Da  59) .  He  found  that  for  a 
deuteron  energy  of  8.2  MeV,  at  reaction  angles  0  =  47°(c.m.)  and  59°(c.m.) 
the  polarization  was  near  -0.10.  By  extrapolating  the  data  below  2  MeV, 
he  predicted  that  for  angles  near  55°(c.m.)  the  polarization  from  2  to 
8  MeV  was  negative,  small  in  amplitude  and  smoothly  varying  with  energy. 

In  the  following  year  Baicker  and  Jones  (Ba  60B)  reported  measurements 
at  a  reaction  angle  of  40° (lab),  which  is  approximately  55°(c.m.)  for  an 
energy  range  of  2  to  4.5  MeV.  Their  results  lie  on  the  curve  predicted 
by  Daehnick  (hereafter  called  the  lower  curve) .  In  the  next  year 
Dubbeldam  and  Walter  (Du  61)  reported  measurements  at  45°(c.m.)  over 
an  energy  range  from  1.9  to  11.0  MeV.  Their  measurements  do  not  agree 
with  the  previous  measurements  and  lie  on  a  curve  which  is  also  a  smooth 
extrapolation  of  the  data  below  2  MeV  (hereafter  called  the  upper  curve) . 
Measurements  have  been  performed  in  other  laboratories  near  this  angle 
and  support  for  either  curve  has  been  reported.  The  measurements  above 
2  MeV  in  Fig.  3-1  are  results  obtained  by  14  different  groups. 

Since  both  Baicker  and  Jones,  and  Dubbeldam  and  Walter  measured 
the  polarization  over  a  range  of  deuteron  energies,  a  closer  investigation 
of  the  difference  of  the  two  sets  of  measurements  was  carried  out  by 
Purser,  Sawers  and  Walter  (Pu  65).  They  used  a  system  similar  to  the  one 
used  by  Dubbeldam  and  Walter  and  measured  the  neutron  polarization  at  a 
reaction  angle  of  45°(c.m.)  over  a  deuteron  energy  range  of  1.9  to  3.7  MeV 
Their  data  are  in  good  agreement  with  the  upper  curve.  Since  Baicker  and 
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Jones  measured  the  neutron  polarization  at  40° (lab),  Purser,  Sawers  and 
Walter  measured  the  polarization  at  this  angle  to  see  if  the  difference 
in  angles  would  account  for  the  discrepancy.  Their  measurements  at  this 
angle  are  slightly  more  negative  than  their  own  measurements  at  45°(c.m.), 
and  statistically  are  still  in  better  agreement  with  the  upper  curve. 

In  order  to  explain  the  discrepancy  between  these  two  sets  of 
measurements,  a  comparison  of  the  experimental  arrangements  employed  by 
both  groups  will  be  discussed.  In  both  cases,  a  deuterium  gas  target 
was  used.  Baicker  and  Jones  used  a  liquid  helium  polarimeter  and  a 
neutron  time-of-flight  spectrometer  set  between  the  target  and  the  neutron 
detector.  The  time  resolution  ranged  from  5  to  10  nano-seconds.  Scintil¬ 
lations  in  the  helium  scattering  volume  were  not  used  to  require  a 
coincidence  with  the  neutron  detector  in  order  to  reduce  the  random 
background.  The  background  reported  in  their  case  ranged  from  40  to  60% 
of  the  total  counts.  Also,  only  one  neutron  detector  was  used,  with  a 
closed  geometry  shielding  consisting  of  a  combination  of  iron,  lead, 
paraffin  and  boric  acid  between  the  target  and  the  polarimeter  and  around 
the  neutron  detector.  Because  of  the  large  bulk  of  the  shielding  materials, 
no  monitor  detector  was  used  and  the  counts  when  the  neutron  detector  was 
placed  in  the  right  and  left  positions  were  normalized  to  the  total  charge 
collected  at  the  target. 

Dubbeldem  and  Walter,  on  the  other  hand,  employed  a  gas  scatterer 
of  54  atmospheres  of  helium  and  5  atmospheres  of  xenon.  The  energy  resolu¬ 
tion  as  obtained  from  the  helium  recoil  scintillation  signals  was  8%  for 
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polonium  a-particles.  A  fast  and  slow  coincidence  was  set  between  the 
polarimeter  and  each  of  the  side  neutron  detectors.  The  background 
reported  in  most  cases  was  not  more  than  20%.  A  solenoid  was  used  to 
change  the  orientation  of  the  neutron  spins  instead  of  moving  the  detectors 
in  the  elimination  of  false  asymmetry.  Open  geometry  shielding  was  employed 
in  this  case.  The  only  shielding  material  in  front  of  the  target  was  a 
brass  collimator  for  the  outgoing  neutrons  aside  from  the  solenoid. 

In  comparing  the  two  systems,  the  most  significant  point  is  that  the 
background  in  the  system  employed  by  Baicker  and  Jones  was  about  three 
times  as  high  as  in  the  system  employed  by  Dubbeldam  and  Walter.  Since 
the  background  was  so  high,  an  over  estimation  of  the  background  in  the 
experiment  of  Baicker  and  Jones  could  result  in  the  magnitude  of  the  neutron 
polarization  being  too  large.  If  this  were  the  case,  the  set  of  measure¬ 
ments  obtained  by  Baicker  and  Jones  would  be  more  negative  than  those  of 
Dubbeldam  and  Walter.  As  the  measurements  of  Baicker  and  Jones  are  more 
negative  than  those  of  Dubbeldam  and  Walter,  their  discrepancy  can  be 
explained  on  account  of  the  high  background  of  the  system  employed  by 
Baicker  and  Jones.  In  the  deuteron  energy  range  of  4  to  9  MeV,  the  polari¬ 
zation  measurements  on  the  lower  and  upper  curves  are  opposite  in  sign,  so 
there  must  be  a  serious  systematic  difference  in  the  systems  employed  by 
the  different  groups  in  this  energy  range  to  obtain  measurements  opposite 
in  sign. 

In  spite  of  the  discrepancy  in  the  neutron  polarization  measurements 
in  the  D(d,n)3He  reaction,  this  reaction  was  chosen  as  a  calibration  of  our 
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system,  because  this  is  the  reaction  on  which  most  polarization  measure¬ 
ments  have  been  performed.  Since  we  also  wanted  to  investigate  the  | 

discrepancy,  the  polarization  of  the  neutrons  from  the  reaction  was 
measured  at  45°(c.m.)  and  at  deuteron  energies  of  3.0,  4.0,  5.0,  and  6.0 
MeV. 

A  solid  target^  (1"  diameter)  was  used  in  this  experiment.  It  was 

a  thin  layer  of  zirconium  (~  500  yg/cm2)  deposited  on  a  tungsten  backing. 

The  deuterium  was  absorbed  in  the  zirconium  layer. 

In  order  to  have  an  unbiased  estimate  of  the  total  number  of  counts 

fi 

under  a  neutron  time-of-f light  peak,  a  computer  program  written  by 

J.  T.  Tepel  (Te  66)  was  used.  Since  a  Gaussian  distribution  is  not 

always  the  correct  shape  in  defining  a  peak  obtained  from  an  experiment, 

this  program  allows  the  input  of  any  desired  reference  peak  shape.  The 

reference  peak  is  defined  by  a  discrete  set  of  (x,y)  co-ordinates  and 

scaling  is  adjusted  by  the  program  to  fit  peaks  (or  overlapping  peaks)  of 

different  heights.  This  program  also  gives  an  assessment  of  the  time 

dependent  background  under  the  peaks  of  the  spectrum.  Linear  and  quadractic 

functions  are  provided  for  background  fitting. 

Our  measurements  of  the  polarization  of  neutrons  from  the  D(d,n)3He 

reaction  at  a  reaction  angle, 0,  of  45°(c.m.)  and  at  deuteron  bombarding 
f 

Isotope  Target  Centre,  Oak  Ridge  National  Laboratory,  Oak  Ridge, 
Tennessee. 

ft 

Adapted  to  run  with  the  IBM  7040  computer  at  the  University  of 
Alberta  by  T.  B.  Grandy  and  D.  A.  Gedcke. 
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energies  of  3.0,  4.0,  5.0,  and  6.0  MeV  are  shown  in  Fig.  3-1,  as  well  as 
many  other  measurements.  They  agree  well  with  the  measurements  by 
Dubbeldam  and  Walter  (Du  61)  and  clearly  support  the  measurements  lying 
on  the  upper  curve.  Our  results  are  also  tabulated  in  Table  3-1  together 
with  the  background  estimation  at  different  deuteron  energies.  The 
increase  in  background  with  increase  in  deuteron  energy  is  due  to  deuteron 
dis association. 


Table  3-1.  Neutron  Polarization  Measurements  from 
the  D(d,n)3He  Reaction  at  9  =  45°(c.m.) 


Ed(MeV) 

<P2  COS(j)> 

Polarization 

Background 

3.0 

.9589 

- (.091+. 026) 

1.3% 

4.0 

.9437 

+(.004+. 020) 

2.7% 

5.0 

.9336 

+(.023+. 021) 

3.5% 

6.0 

.9313 

+  ( . 083+ . 038) 

5.4% 

NOTE:  <P2  cos<J>>  is  the  average  value  of  P2  cost)),  taking 

into  account  the  finite  geometry  of  the  polari- 
meter  and  the  neutron  detectors. 
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CHAPTER  IV 

POLARIZATION  PREDICTIONS  FROM  DIRECT  REACTION  THEORY 

Two  extreme  models  are  usually  employed  to  explain  the  mechanism  of 
a  nuclear  reaction.  The  direct  reaction  mechanism  depicts  the  transition 
from  the  incident  to  the  exit  channel  to  take  place  in  one  step  without 
the  formation  of  an  intermediate  state,  Also  only  a  few  degrees  of 
freedom  are  employed  in  the  description  of  the  colliding  systems.  This 
is  in  contrast  with  the  compound  nuclear  process,  which  assumes  the 
reaction  to  take  place  in  two  steps.  The  incident  particle  is  first 
captured  by  the  target  nucleus  to  form  an  intermediate  compound  nuclear 
state,  which  then  disintegrates  to  yield  the  reaction  products.  The 
distorted  wave  Born  approximation  (DWBA)  has  been  most  commonly  used  in 
the  direct  reaction  calculation  and  many  articles  have  been  written  on 
this  aspect.  The  description  of  the  theory  in  this  chapter  is  based  on 
the  articles  written  by  W.  Tobocman  (To  61) ,  N.  Austern  (Au  63) , 

G.  R,  Satchler  (Sa  64,  Sa  65)  and  R.  H,  Bassel  (Ba  64). 

The  DWBA  theory  of  a  nuclear  reaction  A(a,b)B  is  based  on  the  assump¬ 
tion  that  elastic  scattering  between  a  and  A  in  the  incident  channel  and 
b  and  B  in  the  exit  channel  predominates  over  the  reaction  amplitude. 

Hence  the  reaction  can  be  treated  by  perturbation  theory  and  the  transition 
amplitude  is  given  by 
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In  equation  (4.1),  V  is  the  interaction  potential  responsible  for  the 


—r 

transition.  r  is 
a 


the  displacement  of  a  from  A,  and  r^  the  displacement 
of  b  from  B,  while  J  is  the  Jacobian  of  the  transformation  to  these  relative 


coordinates.  The  functions  x  and  X>-  are  the  distorted  waves,  which  describe 

cL  D 

the  relative  motion  of  the  pair  a,  A(asymptotically  with  relative 

momentum  k^)  before  the  collision,  and  of  the  pair  b^  B  (with  k^)  after  the 

collision.  The  superscript  (+)  or  (-)  denotes  the  outgoing  -  or  incoming  - 

wave  boundary  conditions.  If  the  particles  a  and  b  have  spins  s  and  s  , 

a.  d 

and  the  tatget  and  residual  nuclear  spins  are  J  AND  J  ,  respect ively , we  define 

Ad 

->  -> 

j  ^  JB  "  JA 


s  =  s  s, 

a  b 


and 


->  -* 

^  "  j 


s 


A  partial  wave  expansion  is  used  to  describe  the  distorted  waves  x»  In  the 
presence  of  spin-orbit  coupling,  they  can  be  written  as: 

X^  (k,r)  =  (4n/kr)  )  <(L,s  ,M,m|  J,M-tV)>  <^L,s  ,M+m-m’  ,m'  |  J,M+m^) 

TM+m-m  ’ 


JLM 

M* 


x  i“  xLJ(k,r)Y“  (6k>  *k)YL 


(e  A  ) 

r  r 


(4.2) 


->■ 


The  resultant  angular  momentum,  J  =  L  +  s,  and  its  z-component  (M+m)  are 
conserved  during  scattering  by  the  spin-orbit  potential,  whereas  the 


-> 


individual  z-components  of  L  and  s  are  not.  (Note  that  the  spin  s  used  in 

the  distorted  waves  is  either  s  or  s.  and  should  not  be  confused  with  the 
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spin  s  of  the  captured  particle  as  defined  previously.)  The  terms  for 
which  m  /  m'  allow  for  spin-flips  in  the  reaction.  The  partial  waves, 
Xlj>  are  solutions  of  the  Schrodinger  equation: 


[d2/dr2  +  kz-L(L+l)/r2  -  (2y/h2)U(r)]  x.Lj(k,r)  =  0 


(4.3) 


where  U(r)  represents  the  optical  potentials.  For  large  r,  beyond  the  range 
of  nuclear  potentials,  these  radial  waves  have  the  form: 


XLJ(k,r)  =  |  i  [Hl  (kr)  -  HL(kr) ] exp(iaL) 


(4.4) 


where  =  GL  +  iFL,  and  and  F^  are  the  regular  and  irregular  functions. 

is  the  Coulomb  phase  shift  and  is  the  scattering  matrix  element 

L 

for  the  (L,J)  partial  wave0  The  ru  is  computed  by  numerical  integration 

L 

of  equation  (4.3)  and  then  matching  the  logarithmic  derivative  of  the  wave 
function  to  that  of  equation  (4.4)  at  a  distance, r,  outside  the  nuclear 
field. 

The  optical  model  potentials  are  usually  expressed  as: 


U(r)  =  Vcoul^  "  Vfl(r^  "  iWf2(r)  -  Vgof3(r) 


(4.5) 


Vcoui(r)  as  the  Coulomb  potential.  Assuming  that  a  and  b  are  point  charges 
(ze)  and  A  and  B  are  uniformly  charged  spheres  (z'e),  it  is  given  by: 


Vcoul(r) 
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,  1/3 

where  Rc  is  the  charge  radius  of  the  nucleus  and  is  given  by  r  A 

A  is  the  mass  number  of  the  nucleus, 

V,  W  and  V  are  the  real,  imaginary  and  spin-orbit  well-depth  para¬ 
meters  respectively.  Quite  a  variety  of  shapes  has  been  used  for  the 
radial  functions  fj,  f ? ,  f3  though  the  Saxon-Woods  form  has  been  most 
commonly  used.  They  are: 


f 1 (r)  = 


r-R 


/  -  -R  \ 
1  +  exP  l  ) 


R 


(4.7) 


f2(r)  -  ^ 


4oc 


r-R^ 


1  +  exp("7i) 


(4.8) 


d)  Id 

f  3 (r)  =  -  (— )  -  [f  i  (r)  ]  H  •  s 

J  ^  m  c  J  r  dr  1 

IT 


(4.9) 


where 

and 


RR  “  rR  A 


Rj  =  r^  A 


1/3 

1/3 


Let  us  now  consider  the  nuclear  matrix  element  responsible  for  the 
reaction.  It  can  be  written  as: 


J  <.B.blVlntlA’a>  =  J  <JBMB»abmblVlnt|JAMA>  sama> 
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where  m  =  Mg-MA  +  mb-ma,  The  symbols  bB,  aA  in  the  argument  of  G  denote 
the  dependence  of  G  on  the  various  nuclear  quantum  numbers  other  than  the 
z-components  of  the  spins 0  The  function  G  can  be  defined  by  the  inverted 
form  of  expansion  (4*10) 


G£sJ>m=lSj  ff^I  ..J  <JBMB,  Sb'"blVlnt|jAMA,sama> 


MBMAmbma 


x  (-1) 


sb_mb 


Ja,3  5Ma,Mb-Ma|jb,Mb>  <sa,sb,ma,-mb|s,ma-mb> 


x  <^£,s,m,ma-mb  |  j  ,Mb-Ma^> 


=  A 


Xsj  f£sj  ,m  ^rb  ’  ra^ 


(4.11) 


The  separation  of  G  into  a  spectroscopic  coefficient  A  and  a  form  factor 
f  is  a  matter  of  convenience  in  numerical  computation,  so  that  standard 
form  factors  and  normalization  can  be  used.  For  example,  in  the  study  of 
neutron  hole  states  by  the  (p,d)  or  (3He,a)  reactions,  the  form  factor  is 
the  same  in  both  reactions  but  the  spectroscopic  coefficient,  A  „ ,  is 

X.S  j 

different . 

In  a  (d,n)  stripping  reaction,  a  =  p+n,  b  =  n  and  B  =  A+p.  The 
interaction  potential  can  then  be  written  as: 


V  =  V  -  U 
int  bB  bB 


=  Vnp  +  (VnA  -  unB) 


*  V. 


np 


(4.12) 
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where  UBb  (or  UnB)  is  the  optical  potential  describing  the  scattering  of 
particle  b  (or  n)  by  the  residual  nucleus  B0  It  is  customary  to  take  the 
neutron-proton  interaction  term,  Vnp,  as  the  dominant  term  in  the  inter¬ 
action  potential  and  assume  there  is  a  cancellation  between  the  neutron- 
target  interaction  and  the  optical  potential  in  the  exit  channel,, 

Let  us  now  consider  the  nuclear  matrix  element  of  equation  (4.1) 
once  more.  It  can  be  written  explicitly  as: 


<JBMB,sbmb|Vint|jAMA,Sama>  =  j  %MB  (5aCP  ’  rpA>  *1/2  ^  ^n} 


Vnp  ^n ?p  ’  rnp }  iiu  (  S 5  ^P '  ^ JAMA 


(CA)dCn  d£  d?A 


(4.13) 


where  represents  the  spin  variables  and  the  internal  coordinates  of 
particle  i;  and  tpA  and  rnp  are  respectively  the  distances  between  the 
centres  of  mass  of  the  proton  and  the  target  nucleus,  and  of  the  neutron 
and  the  proton.  In  the  above  expression,  the  interaction,  Vnp,  is  taken 
to  be  scalar,  that  is  only  the  S-state  of  relative  motion  between  the 
neutron  and  the  proton  is  taken  into  consideration,  while  the  D-state  is 
neglected.  Since  the  interaction  is  independent  of  £A,  these  internal 
coordinates  of  the  target  can  be  integrated  over  immedinately .  The  conse¬ 
quent  overlap  with  the  residual  nucleus  B  defines  the  wave  function  <f>  for 
the  captured  proton, 
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<f>0  can  be  expanded  further  as: 
JP 


WV  =  l  <*>t>  VV^WV  y-  (v 


(4.15) 


It  has  been  customary  to  regard  <f>  as  proportional  to  the  shell  model 
single-particle  wave  function  for  the  orbit  (n£j)  so  that 


vv  =  wv 


(4.16) 


where  (r)  is  the  normalized  radial  wave  function  of  the  captured  proton 

and  S  is  the  spectroscopic  factor. 

Putting  equations  (4,14),  (4.15)  and  (4.16)  into  equation  (4.13),  the 
nuclear  matrix  element  can  now  be  written  as: 
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In  the  zero  range  approximation,  one  gets 


V 

np 


(j)  (r  )  =  D06(r  ) 

d  np  u  np 


where 


D0 


(V--a2)  <j>  (r)  dr 


(4,18) 


M  is  the  deuteron  reduced  mass  and  a  is  related  to  the  deuteron  binding 
energy,  e,  by  the  expression 

e  =  iT2or-/M 


The  Hulthen  wave  function  is  usually  chosen  for  <j>^(r). 

Substituting  equations  (4,11)  to  (4,10)  and  comparing  equations  (4.10) 
and  (4,17),  the  spectroscopic  coefficient,  A  and  the  form  factor, 

ioS  J 

f,  .  can  be  written  as: 

Jcsj  ,m 
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(4.19) 


and 


'£sj  ,m 


U„  ,  (r  ) 

£sj  p 


(4.20) 


Note  that  the  Jacobian,  J  in  equation  (4.10)  cancels  out  when  the  zero 
range  approximation  is  carried  out. 

The  transition  amplitude  [equation  (4.1)]  can  now  be  written  as: 

Jimmbma 
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where  m  =  Mg  -  MA  +  mb  -  ma 


The  reduced  amplitude,  6,  is  given  by: 


1/2  llmnbma 
(2j+l)  Ss.  (kb,ka. 


,-V  ->  „ 
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(4.22) 


The  differential  cross  section  of  the  reaction  is  then  given  by: 


dc  =  Vb  XlTl2 

da)  (2Tih2)2  ka  TUa+1.)  (2sa+l) 


(4.23) 


where  y^  and  y^  are  the  reduced  masses  of  the  respective  particles  and  the 
sum  is  over  the  substates  M^,  ma.  Mg  and  mb. 

The  vector  polarization,  Pb ,  of  the  outgoing  particle  is  defined  as: 


P, 
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final 


sb 


^initial^ 


sb 


where  \p .  _  0  -  and  ip-,  ~  are  the  initial  and  final  wave  functions  defining 
initial  final 

the  systems  respectively.  For  convenience  in  numerical  computation,  the 
spectroscopic  coefficient  A  and  the  reduced  amplitude,  3  are  combined  into 
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Using  the  Basel  sign  convention,  the  polarization  of  the  outgoing  particle 


can  now  be  written  as: 


Pb(9) 


1 1  (sb~mb)  (sb+mb+l)  ] J '  "  Im^™^3- 


b£’,iiH-1  ,mb+l  ,m* 
s 'm  J 


£mmbma 

sj 


B£,mmbma*1 
s' j 


(4.24) 


where  the  sums  in  both  the  numerator  and  denominator  are  over  s,  s',  £,£', 
j,  ma,  mb ,  and  m,  In  a  (d,n)  reaction,  Pb(0)  is  the  neutron  polarization 
at  an  emission  angle  of  0. 

In  the  absence  of  spin-orbit  coupling,  the  distorted  waves,  Xmm»  of 
equation  (4.2)  become  diagonal,  since  m=m' .  For  a  deuteron  stripping 
reaction,  equation  (4.24)  can  be  reduced  to: 


where 
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(4.25) 


The  spin  dependent  effect  due  to  the  captured  particle  can  now  be  seen 
from  equation  (4.25),  if  spin-orbit  coupling  in  forming  the  distorted 
waves  is  neglected.  From  equation  (4.25)  one  gets, 
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and 


for  j  =  l  -  j 


(4.26b) 


tt^  is  just  the  expectation  value  of  the  projection  of  the  vector  H  on 
ka  x  lcb,  i.e. 


TT 


l 


and  thus  cannot  exceed  unity.  Since  aj  as  defined  in  equation  (4.25)  is 
also  less  than  unity,  equation  (4.26)  gives  the  limitation  •|Pb(0)|^>‘j  in 


the  absence  of  spin-orbit  coupling.  For  the  case  when  £=0  for  the 
captured  particle,  equation  (4.25)  predicts  the  polarization  Pb(0)  to  be 
zero.  Experimentally,  it  has  been  found  that  the  above  polarization 
predictions  from  DWBA  calculations  without  spin-orbit  coupling  are  not 
true.  Hence  spin-orbit  coupling  must  be  included  in  the  calculation. 

A  DWBA  code  (Ma  64)  written  by  B.E.F.  Macefield  of  Oxford  University, 
England  was  used  in  the  analysis  of  the  neutron  polarization  measurements 
obtained  from  the  28Si(d,n)29P  reaction.  A  discussion  on  this  reaction  is 
given  in  the  next  chapter. 

The  computer  code  was  originally  written  in  Fortran  II  language  and 
used  with  an  IBM  7030  (STRETCH)  computer.  It  computes  the  differential 
cross  section  and  the  polarization  angular  distribution  for  the  outgoing 
particles  in  a  deuteron  stripping  reaction.  Zero  range  and  local 

potentials  are  employed  in  the  calculation;  while  spin-orbit  coupling  in 
the  incident  and  exit  channels ,  as  well  as  for  the  captured  nucleon  is 
included.  The  bound  state  radial  wave  functions  Un^j  (r)  of  the  captured 
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nucleon  are  the  normalized  eigenfunctions  of  the  Schrodinger  equation  with 
a  well  of  the  Saxon-Woods  shape  giving  the  correct  binding  energy.  The 
logarithmic  derivative  of  this  internal  wave  function  is  then  matched 
to  that  obtained  from  the  external  wave  function  (the  Coulomb  function) 
outside  the  nucleus.  The  depth  of  the  Saxon-Woods  well  is  then  adjusted 
successively  until  the  logarithmic  derivatives  are  joined  smoothly 
(Bu  60,  Bu  61) o 

The  shape  of  the  real  part  of  the  optical  potential  is  of  the  Saxon- 
Woods  form.  There  are  three  choices  for  the  shape  of  the  imaginary  part 
of  the  potential.  The  volume  absorption  is  of  the  Saxon-Woods  form,  While  for  sur¬ 
face  absorption  it  can  be  either  Gaussian  or  derivative  of  the  Saxon-Woods 
form.  As  for  the  spin-orbit  potential,  the  shape  is  of  the  Thomas  kind 
(see  equation  4,9).  Also  the  charge  distribution  of  the  nucleus  is  taken 
as  that  of  a  uniformly  charged  sphere. 

The  computer  code  was  translated  into  Fortran  IV  and  used  with  the 
IBM  7040  computer  at  the  University  of  Alberta.  Since  the  memory  core  of 
the  computer  is  only  32  K  words  of  32  bits  each,  the  program  was  divided 
into  three  links.  Also  the  maximum  number  of  partial  waves  used  in  the 
analysis  was  reduced  from  50  to  30. 
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CHAPTER  V 

NEUTRON  POLARIZATION  FROM  THE  28Si(d,n)29P  REACTION 

The  polarization  angular  distributions  of  neutrons  to  the  ground, 

1.38  and  1.96  MeV  states  of  29P  from  the  reaction  z8Si(d,n) z9P  were 
measured  at  a  deuteron  bombarding  energy  of  5.0  MeV.  The  energy  level 
diagram  of  the  nucleus  29P  is  shown  in  Fig.  5-1. 

5 . 1  Experimental  Details 

The  silicon  target  was  made  by  evaporating  a  piece  of  natural  silicon 
onto  a  piece  of  gold  backing  (.005"  thick)  in  vacuum.  An  electron  gun^ 
was  used  in  the  evaporation  process.  The  thickness  of  silicon  in  the 
target  was  about  250  pg/cm2.  An  average  pulsed  deuteron  beam  current  of 
about  2.5  pA  was  used.  The  target  assembly  was  cooled  by  a  rapid  flow 
of  air. 

In  order  to  eliminate  false  asymmetry,  the  two  neutron  detectors 
were  interchanged  in  their  right-left  positions  for  the  measurement  at 
every  angle.  Typically  it  took  3-1/2  hours  to  accumulate  400  counts  under 
the  smallest  time-of-flight  peak  of  interest.  Hence  for  the  measurement 
at  one  angle,  it  took  about  7  hours  of  counting  time. 

Fig.  5-2  shows  a  neutron  time-of-flight  spectrum  of  the  28Si(d,n)29P 
reaction  at  a  neutron  emission  angle  of  60°  (lab) .  The  neutron  flight 
path  was  2  metres.  A  scale  of  the  neutron  flight  time  is  shown  as  a 

^Model  CP-500,  Andar  Corporation,  Mountain  View,  California 
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Fig,  5-1 


The  energy  level  diagram  of  the  nucleus 
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Fig.  5-2 


A  neutron  time-of-f light  spectrum  of  the  Si(d,n)  P 
reaction  with  E^=5.0  MeV,  0=60° (lab)  and  flight  path  =2m. 
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vertical  axis  on  the  right-hand  side  of  the  figure.  The  two  fluorine 
states  in  the  time-of-flight  spectrum  were  due  to  contamination 

on  the  surface  of  the  silicon  target.  Carbon  contamination  was  found 
to  be  very  small. 

The  neutron  time-of-flight  spectra  of  this  reaction  as  recorded 

by  the  polarization  system  have  been  shown  in  Fig.  2-13  of  chapter  2. 

They  correspond  to  the  true  and  random  time-of-flight  spectra  of 

neutrons  scattered  to  the  right  or  left  by  the  polarimeter.  The 
f 

computer  code  written  by  Tepel  (Te  66)  was  used  to  find  the  total 
number  of  counts  under  the  different  peaks. 

5. 2  Ground  State  of  29P 

Fig.  5-3  shows  the  angular  distribution  and  the  polarization  angular 
distribution  of  neutrons  to  the  ground  state  of  29P.  The  angular  distri¬ 
bution  was  measured  by  W.  G.  Davies  j2t  al_  (Da66  A)  at  the  University  of 
Alberta.  The  curves  in  the  two  graphs  are  drawn  through  the  experimental 
points  to  show  the  trend  of  the  distributions. 

Since  £p  =  0  for  the  capture  of  the  proton  by  the  28Si  nucleus  to 
form  the  ground  state  of  29P,  the  polarization  angular  distribution  of 
neutrons  corresponding  to  this  state  is  expected  to  obey  the  derivative 
rule  (Ro  59,  Bi  60,  Jo  62).  In  equation  form,  this  rule  states: 

P  ( 6)  « 

*j* 

See  end  of  chapter  3  for  a  brief  description  of  the  computer  code. 


ae  a(e). 


/o(e) 


for  £  =  0 


(5.1) 
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Fig„  5-3  Angular  distribution  and  polarization  angular  distribution 

29 

of  neutrons  to  the  ground  state  of  P 


A 

' 

— 
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The  behaviour  of  the  polarization  angular  distribution  from  our 
measurements  as  compared  to  the  derivative  rule  can  be  seen  from  Fig*  5-3. 
It  increases  negatively  from  zero  in  accordance  with  the  sign  of  the  slope 
of  the  angular  distribution.  It  crosses  zero  first  at  35°  which  is  the 
minimum  of  the  angular  distribution.  It  then  increases  positively  and 
is  again  in  agreement  with  the  sign  of  the  slope  of  the  angular  distri¬ 
bution.  It  crosses  zero  again  at  72°  which  is  about  8°  from  the  secondary 
maximum  of  the  angular  distribution.  It  finally  crosses  zero  at  110°  which 
is  the  second  minimum  of  the  angular  distribution.  As  for  the  magnitude 
of  the  polarization,  it  is  also  clear  that  it  gets  larger  at  the  turning 
points  as  the  angle  increases.  Hence,  qualitatively  at  least,  the 

distribution  from  29P  ,  agrees  well  with  the  derivative  rule  as  defined 

gnd 

by  equation  (5-1). 

According  to  equation  (4-25)  in  chapter  4,  the  polarization  for  £  =  0 
should  be  equal  to  zero,  if  spin-orbit  coupling  is  neglected.  This  is 
surely  not  true  in  the  29P  ground  state  measurements.  This  indicates  that 
spin-orbit  coupling  must  be  included  in  the  calculation  of  polarization. 

A  tablulation  of  the  results  corresponding  to  the  ground  state  of 
29P  is  shown  in  Table  5-1. 

5. 3  The  1.38  and  1.96  MeV  States  of  29P 

The  angular  distributions  and  the  polarization  angular  distributions 
of  neutrons  corresponding  to  the  1.38  and  1.96  MeV  states  are  very  useful 
in  the  study  of  j-dependence  of  the  captured  nucleon.  The  £-value  for  the 
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■ 


. 
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Table  5-1 


o  Polarization  of  neutrons  from  the 
28Si(d,n)29P  ^  reaction  at  = 
5.0  MeV 


0 (lab) 
deg. 

0  (  c .  m . ) 
deg. 

<P2COS(j)> 

Polarization 

10 

10.47 

.9519 

-(.010  +  .017) 

15 

15.70 

.9563 

-(.038  +  .024) 

20 

20.93 

.9583 

-(.028  +  .022) 

30 

31.36 

.9605 

-(.139  +  .037) 

40 

41.75 

.9619 

+(.204  +  .022) 

50 

52.08 

.9631 

+(.133  +  .018) 

60 

62.36 

.9641 

+(.130  +  .051) 

70 

72.56 

.9651 

+(.002  +  .023) 

80 

82.65 

.9659 

-(.163  +  .028) 

90 

92.72 

.  9666 

-(.178  +  .040) 

100 

102.68 

.9672 

-(.190  +  .047) 

110 

112.56 

.9677 

+(.045  +  .034) 

120 

122.36 

.9681 

+(.209  +  .029) 

' 


ex 


■ 
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captured  proton  is  2  in  these  two  states.  For  the  1.38  MeV  state 

13  15 

j  =  £  -  =  ~2  >  while  for  the  1.96  MeV  state  j  =  l  +  —  =  —  . 

Lee  and  Schiffer  (Le  64A)  noticed  in  their  survey  of  (d,p)  stripping 

reactions  that  for  £n  =  1  and  j  =  y  states  the  angular  distribution 

shows  a  minimum  between  100°  and  135°  (c.m.),  while  that  for  £n  =  1  and 
3“ 

j  =  —  does  not  show  a  minimum  in  this  angular  region.  They  attribute 
this  difference  to  the  j -dependence  of  the  captured  neutron,  and  suggest 
that  such  an  indication  in  the  angular  distribution  might  be  useful  in 
assigning  j-values. 

On  the  basis  of  DWBA  calculations,  the  spin-orbit  interaction  of  the 
captured  nucleon  is  attractive  for  the  j  =  £  +  y  transition,  while  the 
opposite  effect  occurs  in  the  j  =  £  -  y  transition.  However,  such  effects 
have  been  estimated  to  be  too  small  to  change  the  shape  of  the  angular 
distribution,  although  the  magnitude  of  the  cross  section  is  affected  (Le  64B) . 

The  j -dependence  for  the  1.38  and  1.96  MeV  states  of  29P  is  shown  in 
Fig.  5-4.  The  4.0,  5.0  and  6.0  MeV  angular  distributions  were  measured  by 
W.  G.  Davies  et_  al  (Da  66B),  while  the  7.0  and  9.0  MeV  angular  distributions 
were  measured  by  S.  G.  Buccino  et_  a_l  (Bu  66).  The  curves  are  drawn  through 
the  experimental  points  to  show  the  trends  of  the  distributions  and  they 
are  normalized  at  the  peak  to  show  the  j -dependence  of  the  two  states.  For 
angles  beyond  the  main  stripping  peak,  the  angular  distribution  for  the 
1.38  MeV  state  (j  =  £  -  y  )  drops  faster  than  that  for  the  1.96  MeV  state 
(j  =  £  +  y  ), except  at  =  4  MeV,  where  compound  nucleus  effects  may  play 
a  major  role.  This  difference  is  observed  to  be  more  prominent  at  higher 


deuteron  energies. 
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Fig,  5 


4  The  j-dependence  effect  as  seen  from  the  angular  distributions 


of  the  1.38  and  1.96  MeV  states  of 
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The  polarization  measurements  corresponding  to  the  1.38  and  1.96  MeV 
states  are  shown  in  Fig.  5-5.  The  angular  distribution  is  shown  on  the 
top  of  the  figure  for  comparison.  Again,  the  curves  are  drawn  through 
the  experimental  points  to  show  the  trends  of  the  distributions.  At  such 
a  deuteron  bombarding  energy  (5.0  MeV),  the  difference  in  the  angular 
distributions  of  these  two  states  is  small,  but  the  difference  in  the 
polarization  distributions  is  quite  large,  especially  at  forward  and 
backward  angles.  Hence  polarization  is  a  more  sensitive  way  of  showing 
this  j-dependent  effect. 

According  to  the  sign  rule  as  predicted  by  Newns  (Ne  53) ,  the 

polarization  should  be  positive  around  the  main  stripping  peak  for  j  =  £  + 

-j  and  negative  for  j  =  £  -  -j  transitions.  Very  often  this  rule  is  not 
obeyed.  Instead,  D.  W.  Miller  (Mi  65)  noticed  that  for  angles  just  beyond 
the  main  stripping  peak,  the  polarization  is  small  and  fluctuates  around 
zero  for  j  =  £  +  y  and  is  negative  and  larger  for  j  =  £  -  —  from  his 
survey  of  (d,p)  polarization  measurements.  In  our  polarization  measure¬ 
ments,  the  polarization  for  the  1.96  MeV  state  (j  =  £  +  -  )  is  small  and 

mostly  negative  for  the  angular  range  of  40°  to  70°,  while  that  for  the 

1.38  MeV  state  (j  =  £  -  ~)  is  negative  and  definitely  larger.  Hence 
this  observed  trend  in  (d,p)  polarization  measurements  applies  to  these 
two  levels  as  well. 

The  negative  dip  near  15°  for  the  1.96  MeV  state  is  also  observed 
in  our  polarization  measurements  of  the  x 60(d,n) 1 reaction.  The 

5+  1 

j-value  for  the  captured  proton  in  this  case  is  also  —  ( j  =  £  +  ) • 
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Fig.  5 


5  Angular  distributions  and  polarization  angular  distributions 


of  neutrons  to  the  1.38  and  1.96  MeV  states  of 


29 
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Hence  it  would  be  interesting  to  see  if  other  measurements  for  j  = 
transitions  would  follow  this  trends 

Tabulations  of  the  polarization  measurements  for  the  1.38  and  1.96 
MeV  states  are  shown  in  Tables  (5-2)  and  (5-3)  respectively. 


Table  5-2.  Polarization  of  neutrons  from  the 
28Si(d,n) 29P],  30  reaction  at  Ej  = 
5. 0  MeV 


6 (lab) 
deg. 

0 (c.m. ) 
deg. 

(Pjcos  4>^> 

Polarization 

10 

10.55 

.9526 

+(.024  +  .032) 

15 

15.82 

.9567 

-(.175  +  .041) 

20 

21.09 

.9583 

-(.028  +  .022) 

30 

31.59 

.9595 

-(.083  +  .019) 

40 

42.04 

.9595 

-(.099  +  .013) 

50 

52.43 

.9589 

-(.136  +  .022) 

60 

62.75 

.9579 

-(.146  +  .032) 

70 

72.99 

.9564 

-(.185  +  .034) 

80 

83.13 

.9545 

-(.253  +  .042) 

90 

93.18 

.9523 

-(.153  +  .027) 

100 

103.13 

.9499 

-(.133  +  .035) 

110 

112.99 

.9473 

-(.146  +  .029) 

120 

122.75 

.9466 

-(.186  +  .023) 

LO|  CM 


' 


- 


Table  5-3.  Polarization  of  neutrons  from  the 
28Si(d,n) ^3Pi . gg  reaction  at  = 
5.0  MeV 


0 (lab ) 
deg. 

0 (c.m. ) 
deg. 

(P  2cos  <J>^ 

Polarization 

10 

10.60 

.9350 

-(.155  +  .058) 

15 

15.89 

.9388 

-(.296  +  .070) 

20 

21.18 

,9401 

-(.249  +  .051) 

30 

31.  73 

.9403 

-(.153  +  ,048) 

40 

42.22 

.9390 

-(.043  +  .022) 

50 

52.64 

.9366 

-(.064  +  .035) 

60 

62.99 

.9335 

+(.045  +  .028) 

70 

73.24 

.9295 

-(.075  +  .037) 

80 

83.40 

.9250 

-(.030  +  .031) 

90 

93.45 

.9198 

-(.111  +  .040) 

100 

103.40 

.9143 

-(.280  +  .060) 

110 

113.24 

.9085 

-(.381  +  .036) 

120 

122.99 

,9028 

-(.410  +  .056) 

' 
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5-4  DWBA  Predictions 

A  comparison  was  made  between  the  DWBA  code  written  by  Macefield 

(Ma  64)  and  the  code  "JULIE"  written  by  Bassel,  Drisko  and  Satchler  (Ju  66). 

40  41 

The  output  data  for  the  Ca  (d,p)  Ca  ^  reactaon  at  E^=12  Mev  from  the  two 
codes  are  shown  in  Fig  (5-6).  The  shape  of  the  angular  distributions 
compares  very  well.  "JULIE"  estimates  the  cross  section  at  the  main 
stripping  ,peak  to  be  about  10%  higher  than  that  estimated  by  Macefield. 

The  polarization  outputs  do  not  compare  as  well.  However,  with  the  present 
statistical  accuracy  in  polarization  measurements,  one  curve  is  as  useful 
as  the  other  in  comparing  theory  with  experiment  for  the  present,  though 
of  course  the  question  of  accuracy  of  the  Macefield  program  must  be 
investigated  more  thoroughly  when  the  new  University  of  Alberta  computer 
is  operational. 

The  code  written  by  Macefield  has  been  used  to  calculate  the 

polarization  of  neutrons  corresponding  to  the  ground,  1.38  and  1.96  MeV 
29  28  29 

states  of  P  emitted  in  the  Si(d,n)  P  reaction.  The  deuteron  optical 


parameters  are  taken  from  the  results  of  Perey  and  Perey  (Pe  66) .  They 

28 

analysed  the  elastic  scattering  angular  distributions  of  deuterons  on  Si 
over  an  energy  range  of  10.0  to  11.8  MeV.  They  noticed  that  the  inclusion 
of  a  spin-orbit  potential  improved  the  fits  at  large  angles.  From  the 
automatic  search  code,  they  arrived  at  two  sets  of  opitcal  parameters,  one 
with  the  real  well  depth  at  65.95  MeV  and  the  other  at  188.5  MeV.  These 
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Fig.  5-6  A  comparison  between  "JULIE"  and  the  DWBA  code  written  by 
Macefield 
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two  sets  were  used  in  an  attempt  to  fit  our  measurements,  and  a  better 
fit  is  found  with  the  second  set.  As  for  the  neutron  optical  parameters, 
the  set  of  optical  parameters  found  by  Rosen  (Ro  65)  for  the  energy  range 
of  0.2  to  24  MeV  is  used.  Buck  and  Perey  ^ also  suggested  a  set  of 
neutron  optical  parameters,  but  this  set  is  found  to  give  a  fit  not 

as  good  as  the  set  by  Rosen.  As  for  the  captured  proton,  the  set  of 

optical  parameters  as  used  by  Davies  (Da  66A)  has  been  used  in  the 

calculation  of  the  ground  state.  If  the  same  set  of  parameters  is  used 

for  the  1.38  and  1.96  MeV  states,  the  real  well  depth,  which  is  auto¬ 
matically  adjusted  by  the  DWBA  code,  is  more  than  100  MeV  for  these  two 
states.  For  this  reason,  the  real  radial  parameter,  r  ,  and  the  real 
diffuseness,  a  ,  are  made  larger  to  keep  the  real  well  depth  below  100  MeV. 

A  tabulation  of  the  optical  parameters  for  the  deuteron,  neutron  and 
captured  proton  is  shown  in  Table  (5-4) .  The  form  of  the  optical  potentials 
are  defined  by  equations  (4-5)  to  (4-9)  in  Chapter  4. 

The  DWBA  predictions  for  the  angular  distributions  and  polarization 

angular  distributions  of  neutrons  corresponding  to  the  ground,  1.38  and 

29 

1.96  MeV  states  of  P  are  shown  in  Figs. (5-7),  (5-8)  and  (5-9)  respectively. 

The  predictions  for  the  angular  distributions  give  the  right  shape  near 
the  main  stripping  peak  and  they  are  too  low  for  angles  beyond  the  main  stripping 
peak.  The  polarization  prediction  for  the  ground  state  is  fairly  good  at 
forward  angles,  but  at  the  backward  angles,  it  becomes  too  large.  For 


^reported  in  reference  Ro  65. 
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O  O  OQ 

Table  (5-4)  OPTICAL  MODEL  PAPAMETERS  USED  IN  THE  Si(d,n)  REACTION 

AT  DEUTERON  ENERGY  OF  5.0  MeV 


Deuteron 

Neutron 

V 

188.05 

47.07 

w 

16.06 

5.75 

vso 

11.21 

5.05 

rR 

.651 

1.25 

ri 

1.518 

1.25 

rcoul 

1.3 

— 

aR 

1.047 

.  65 

ai 

0.564 

.  7 

CAPTURED 

PROTON 

29p  , 

gnd 

29pl. 38 

29pl. 96 

V 

51.8 

73.0 

67.3 

V 

so 

— 

5.0 

5.0 

r 

R 

1.25 

1.5 

1.5 

rcoul 

1.25 

1.25 

1.25 

aR 

.5 

1.0 

1.0 

NOTE:  All  potentials  are  in  MeV  and  all  lengths  are  in  fermis 
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Fig.  5-7  DWBA  predictions  for  the  angular  distribution  and  polarization 

angular  distribution  of  neutrons  corresponding  to  the  ground 
29 

state  of  P 


POLARIZATION 
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Fig  5-8  DWBA  predictions  for  the  angular  distribution  and  polarization 

angular  distribution  of  neutrons  corresponding  to  the  1.38  MeV 
29 

state  of  P 


100 
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Figo  5-9  DWBA  predictions  for  the  angular  distribution  and  polarization 

angular  distribution  of  neutrons  corresponding  to  the  1„96  MeV 
29 

state  of  P 


POLARIZATION 


CP 
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thel°38  MeV  state,  the  DWBA  prediction  does  not  agree  with  the  polarization 
measurements  at  all.  For  the  1096  MeV  state,  it  gives  the  features  of  the 
polarization  measurements  up  to  about  80°,  and  beyond  that  it  even  becomes 
opposite  in  sign. 

It  should  be  emphasized  that  there  are  no  adjustments  in  the  optical 
parameters  of  the  incident  and  exit  channels  to  improve  the  fits.  They 
are  just  the  optical  parameters  as  determined  from  elastic  scattering 
experiments.  It  is  quite  unfortunate  that  there  are  no  optical  parameters 
from  elastic  scattering  data  near  a  deuteron  energy  of  5  MeV  and  we  have 
had  to  use  the  parameters  of  Perey  and  Perey  for  higher  deuteron  energy. 

5-5  Compound  and  Direct  Effects 

The  yield  functions  of  neutrons  corresponding  to  the  ground,  1.38 

29  28  29 

and  1.96  MeV  states  of  P  from  the  reaction  Si(d,n)  P  were  measured  at 
20°,  60°  and  120°  (lab)  over  the  deuteron  energy  range  of  4,64  to  5  =  50  MeV. 
They  are  shown  in  Figs, (5-10),  (5-11)  and  (5-12).  All  the  yield  functions 
fluctuate  quite  a  bit.  They  are  different  for  different  states,  and  even 
for  the  same  state  they  are  different  at  different  angles.  One  obvious 
conclusion  that  can  be  drawn  from  the  fluctuations  in  the  yield  is  that 
apart  from  the  direct  reaction  contribution,  a  compound  nuclear  effect  is 
present  in  this  energy  region.  This  might  explain  the  poor  agreement 
between  the  DWBA  predictions  and  the  experimental  measurements  discussed 
in  the  last  section.  Another  point  that  is  worth  mentioning  is  that  the 

28 

t  Elastic  scattering  analysis  of  5.0  MeV  deuteron  on  Si  is  now  in  process 
by  D.  Gurd  of  this  laboratory. 
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Fig.  5-10  Yield  functions  of  neutrons  corresponding  to  the  ground 

29 

state  of  P 
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Figc  5-11  Yield  functions  of  neutrons  corresponding  to  the  1=38  MeV 


29 

state  of  P 
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Fig„  5-12  Yield  functions  of  neutrons  corresponding  to  the  1.96  MeV 

29 

states  of  P 
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the  measured  polarization  angular  distribution  corresponding  to  the  ground 
29 

state  of  P  is  in  good  qualitative  agreement  with  the  derivative  rule. 

Since  this  rule  is  derived  by  means  of  an  approximate  DWBA  calculation,  this 
agreement  may  indicate  that  the  direct  reaction  contribution  is  still 
predominant  in  this  energy  ..region.  It  is  difficult  to  estimate  the 
contribution  from  compound  nuclear  processes.  A  method  in  the  separation 
of  the  direct  reaction  contribution  and  the  compound  nuclear  effect 
is  described  by  Austern  in  Chapter  V.D.  of  Fast  Neutron  Physics,  Part  II 
(Au  63) . 
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APPENDIX  I 

CALCULATION  OF  THE  n-4He  ANALYSING  POWER 

In  this  section  a  description  of  a  computer  program  written  for  the 
calculation  of  the  n-4He  analysing  power  is  given.  In  the  calculation, 
a  point  target  was  assumed  while  the  finite  geometry  of  the  scattering 
volume  of  the  polarimeter  and  the  neutron  detector  was  taken  into  con¬ 
sideration,  Since  the  maximum  angle  subtended  by  the  scattering  volume 
at  a  distance  of  2  metres  from  the  target  was  only  +  1,5°,  the  variation 
of  cross  section  and  polarization  of  the  nuclear  reaction  over  this 
angular  range  was  not  taken  into  consideration.  The  program  was  written 
in  Fortran  IV  language  and  used  with  the  IBM  7040  computer  at  the  University 
of  Alberta. 

Recently  a  computer  program  has  been  written  by  G.  M,  Stinson  and 
S.  M.  Tang  of  this  laboratory  to  include  the  double  scattering  effect  in 
the  calculation  of  the  analysing  power.  With  the  inclusion  of  this  effect 
the  analysing  power  is  reduced  by  about  10  %.  Since  the  difference  is 
within  the  statistical  uncertainty  of  the  polarization  measurements,  the 
analysing  power  used  in  this  thesis  is  calculated  by  considering  single 
scattering  alone. 


■ 


. 


94 


A1 . 1  Principle  of  Calculation 

Fig,  Al-1  shows  a  schematic  drawing  of  the  polarization  experiment 
configuration. 


Fig.  Al-1.  A  schematic  drawing  of  the 
polarization  configuration. 

In  this  configuration,  the  target,  the  centre  of  the  scattering 
volume  and  the  centres  of  the  two  neutron  detectors  are  in  a  horizontal 
plane.  Let  the  distance  between  the  target  and  the  centre  of  the  scatter¬ 
ing  volume  be  S]_q  and  the  distance  between  the  centre  of  the  scattering 
volume  and  the  centre  of  the  front  face  of  the  neutron  detector  be  S20* 

The  scattering  volume  is  placed  at  an  angle  of  0]_q  with  respect  to  the 
incident  deuteron  beam,  which  is  assumed  to  be  horizontal,  while  the 
neutron  detectors  are  placed  at  an  angle  of  +  ©20  with  respect  to  the 
line  joining  the  target  and  the  centre  of  the  scattering  volume.  For 

the  purpose  of  convenience  in  numerical  computation,  the  following  co- 

i 

ordinate  systems  are  introduced. 

^The  right  hand  rule  is  used  in  defining  all  the  co-ordinate  systems 
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(I)  Deuteron  reaction  frame  of  reference  with  origin  at  the  target . 

Z±  axis  -  along  the  incident  deuteron  direction  and  horizontal 
Xi  axis  -  90°  from  Z\  axis  and  horizontal 
Yi  axis  -  normal  to  plane,  i.e.  vertical 

(II)  This  frame  of  reference  is  obtained  by  rotating  about  the  axis 
through  an  angle  of  0io>  such  that 

Z2  axis  -  along  the  line  joining  the  target  and  the  centre  of  the 
scattering  volume 

X2  axis  -  90°  from  Z2  axis  and  horizontal 
Y2  axis  -  same  as  Yj  axis  and  vertical 

(III)  Neutron  scattering  frame  of  reference 

This  frame  of  reference  is  obtained  by  translating  along  Z2  axis 
through  a  distance  of  S^q,  such  that 

Z3  axis  -  along  the  line  joining  the  target  and  the  centre  of  the 
scattering  volume 

X3  axis  -  90°  from  Z3  axis  and  horizontal 
Y3  axis  -  normal  to  Z3X3  plane  and  vertical 

(IV)  This  frame  of  reference  is  obtained  by  rotating  about  Y3  axis 
through  an  angle  o  f  ®20  >  such  that 

Z4  axis  -  along  the  line  joining  the  centres  of  the  scattering 
volume  and  the  neutron  detector 

X4  axis  -  90°  from  Z4  axis  and  horizontal 
Y4  axis  -  normal  to  Z4X4  plane  and  vertical 
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(V)  This  frame  of  reference  is  obtained  by  translating  along  Z4  axis 
through  a  distance  of  S20J  such  that 

Z5  axis  -  along  the  line  joining  the  centres  of  the  scattering 
volume  and  the  neutron  detector 

X5  axis  -  90°  from  Z5  and  horizontal 

Y5  axis  -  normal  to  Z5Y5  plane  and  vertical 

For  the  finite  geometry  calculation  of  the  analysing  power  of  n-^He, 
both  the  scattering  volume  and  the  neutron  detectors  are  divided  into 
small  rectangular  cells,  and  the  scattering  or  the  detection  of  the 
neutrons  are  assumed  to  occur  at  the  centres  of  the  cells,,  Since  the 
scattering  volume  is  a  cylinder  of  2  inches  in  diameter  and  4  inches  in 
height,  the  boundary  conditions  for  the  co-cordinates ,  (X3,  y3,  Z3) ,  of 
the  centres  of  each  of  the  cells  are: 


-(2 


and 


where  AY 3 


AY  3  AY  3 

x  2.54  -  ——  )  <_  y  3  (2  x  2.54  -  — —  )  cm 


z^  +  x§  <_  (2«54)2  cm2 


is  the  cell  dimension  in  the  Y3~axis  direction. 


As  for  the  neutron  detector,  which  is  3-1/2  inches  in  diameter  and  2  inches 
in  length,  the  boundary  conditions  for  the  co-ordinates  (X5,  y5,  Z5)  of  the 
centre  of  each  of  the  cells  are: 


AZ5 

<  z5  <_  (2  x  2.54 


AZ5 
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x5  +  y  5  C^T-  x  2. 54)  2  cm2 

where  AZ5  is  the  cell  dimension  in  the  Z5-axis  direction* 

Let  us  now  follow  one  of  the  neutron  events  in  the  system  and  calculate 
the  different  angles  required  in  the  numerical  computation*  The  direction 
cosines  of  the  incident  deuteron  direction  with  respect  to  System  II  are: 

Aj  :  Bj  :  C]_  =  -  sinS^o  :  0  :  cosQ^q 

Since  (X3,  y3,  Z3)  are  the  co-ordinates  of  the  centre  of  a  cell  in  the 
scattering  volume  with  respect  to  System  III,  then  with  respect  to  System 
II  the  co-ordinates  are: 

(x3 >  y 3 >  S10  +  23) 

The  direction  cosines  of  the  line  joining  the  target  and  the  centre  of  the 
cell  in  System  II  are: 

x3  y3  s10+z3 

a2  :  B2  :  C3  =  —  :  —  :  —5— 

where  r  2  ,  2  .  ,c  .  \2i1/2 

G  =  [x|  +  y|  +  (S10  +  z3)z] 

These  are  just  the  direction  cosines  of  the  outgoing  neutron*  Hence  the 
angle,  0,  between  the  incident  deuteron  beam  and  the  outgoing  neutron  is 
given  by : 


cos  0 ]_  =  A1A2  ®l^2  ^1^2 


(Al-1) 
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The  direction  cosines  of  the  normal  to  the  reaction  plane  are: 
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B  ^2  -  B  ?C 


B 


2W 


C,  = 


D 


A2C1  ~  A1C2 
D 


Al®2  —  A9B 


2D  1 


D 


(Al-2) 


where 


D  =  [ (B x C 2  ~  B2C! ) ^  +  (A2C !  -  A!C2  )2  +  (A^  -  A2B!)2] 


1/2 


With  respect  to  System  V,  the  co-ordinates  of  the  centre  of  a  cell  in  the 
neutron  detector  are  (X5,  y5,  Z5)  With  respect  to  System  IV  the  co-ordinates 
are: 

(x5 »  y5>  S20+  z  5) 

With  respect  to  System  III,  the  co-ordinates  are: 

x 5  COS  e2o+  (S20  +  z 5) sin  e2o 

y  5 

and  -X5  sin  @20+  (S20  +  Z5)cos  ©20 

With  respect  to  System  II,  the  co-ordinates  are: 


=  X5  cos  020"*"  (S20  +  Z5)sin  @20 

H2  =  ys 

H3  =  -x5  sin  02O  +  (S20  +  z 5) cos  02o+  S10 


where  (H]_,  H2 ,  H3)  are  the  co-ordinates  of  the  centre  of  the  cell  in  the 
neutron  detector  with  respect  to  System  II. 

The  direction  cosines  of  the  straight  line  joining  the  centres  of  the  two 
cells  i.e.(x3,  y 3,  S^q  +  z3)  in  the  scattering  volume  and  (Hj,  H2,  H3)  in 
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the  neutron  detector  are: 


A4  :  B4 


h2  -  y  3 
E 


h3  “  Cs  1 0  +  23) 
E 


where  E  +  {(Hi  -  x3)2  +  (H2  -  y3)2  +  [H3  -  (S 1 0  +  z3)]2}1//2 

These  also  are  the  direction  cosines  of  the  scattered  neutron  direction. 
Hence  the  angle,  02,  between  the  incident  neutron  and  the  scattered 
neutron  is  given  by : 


COS  02  =  A2A;4  +  B2B4  +  C2C4 


(Al-3) 


The  direction  cosines  of  the  normal  to  the  scattering  plane  are: 


Ac  :  Be  :  C  c  = 


B2C4  —  B4C;?  A4.C2  —  A2C4  A2B4  —  A4B2 


(Al-4) 


where 


l/2 

F  =  [  (B2C4  —  B4C2) 2  +  (A4C2  —  A2C4)2+  (A2B4  —  A4B2)Z-  ] 


From  equations  (Al-2)  and  (Al-4),  the  angle,  <J>,  between  the  normals  of  the 
reaction  and  the  scattering  planes  is  given  by 

cos  <j)  =  A3A5  +  B 3B 5  +  C3C5  (Al-5) 

Because  of  the  finite  size  of  the  scattering  volume  and  the  neutron  detectors, 
cos  cp  is  equal  to  1  only  for  the  case  when  the  production  and  the  scattering 
of  the  neutrons  occur  in  the  same  plane.  In  general,  if  the  angle  between 
n3  and  n2  is  <f>  (see  equation  1-1  in  Chapter  1)  ,  the  right-left  asymmetry,  e, 
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is  given  by 

e  =  P1P2  cos  <j> 

instead  of  the  expression  in  equation  (1-2).  Since  the  differential  cross 
section,  a,  of  the  scattered  neutron  is  given  by 

a  =  |  g  |  2  +  |  h  |  2  , 

equation  (1-3)  can  be  written  as: 

P2a  =  -2Im(g*h) 

In  the  computer  program,  P2C  cos  <p  and  a  are  calculated  for  every  combina¬ 
tion  of  the  cells  in  the  scattering  volume  and  the  neutron  detector.  The 
average  value  of  P2  cos  <f>  is  then  given  by: 


N 

l  (P2a  cos  (p)± 

i=l _  (Al-6) 

N 

l  °± 

i=l 

where  N  is  the  total  number  of  the  possible  combinations  of  the  cells.  In 
this  way  the  average  analysing  power,  ^2  cos  is  calculated. 

For  the  completeness  in  describing  the  computation  process  in  the 
program,  some  algebraic  expressions  will  be  given  in  the  following  para¬ 
graphs  . 

When  6^  =  0  for  i  >  2,  g  and  h  (see  equation  (1-3))  can  be  written  as 


<^P2  cos  4> 


' 

b  ■  .  1  if  o  ,  liar .♦*  >i  : 

>n  v  g  •"  ■  • 


>v 
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X  i&0  4-  10 1 

g  =  T"  (  P o  (cos0)  [sin6o  e  ]  +  Pj  (cos0)  [2sin6  i  e  +  sin6  j  e  ] 


16 


+ 


16' 


_l_  i6  2  _  i62 

+  P2(cos0) [3sin62  e  +  2sin62  e  ]} 


15 1 


i(«t  +'s7> 


h  =  —  { Pq (cos0)sin6o  e  +  P J (cos 0) sin(6 j  -  6j)  e 


+ 


+  P2(cos0)sin(6  -  6„)  e 


i(62  +  62) 


} 


Separating  g  into  its  real,  GR,  and  imaginary,  GI,  parts. 


GR  =  (GRO)  +  (GR1) cos0  +  (GR2)  cos20 

1  3  +  +  -  - 

where  GRO  =  —  {sin60  cos6q  -  —  sin62  cos62  -  sin62  cos 6 2  } 


GR1  =  ^  {2sin6^"  cos6+  +  sin  6^  cosd^^  } 


1  9  +  +  _  _ 

GR2  =  — {  —  sin62  cos60  +  3sin62  cos62  } 


and  GI  =  (GIO)  +  (Gil)  cos0  +  (GI2)  cos20 

where  GIO  =  7-  {sin260  -  7  sin26+  -  sin269  } 
k  u  2  2  ^ 


Gil  =  f  {2sin26t  +  sin26,  } 

k  1  1 


GI2  =  £  {-|  sin262  +  3  sin262  } 


Similarly,  separating  h  into  its  real,  HR,  and  imaginary,  HI,  parts, 


HR  =  (HR1)  sin0  +  (HR2)  sin  20 
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where 


and 


HI  =  (HI1)  sine  +  (HI2)  sin28 


where 


These  phase  shifts  have  been  expressed  as  power  series  in  neutron  energies 
in  Section  1-4  of  Chapter  I.  A  subroutine  has  been  written  to  calculate 
the  energies  of  the  neutrons  emitted  at  different  angles  from  a  nuclear 
reaction.  This  subroutine  also  calculates  the  centre-of-mass  angle  corres¬ 
ponding  to  the  angle,  02  (equation  (Al-3)),  between  the  incident  neutrons 
and  the  scattered  neutrons,  since  the  angles  in  the  expressions  for  g  and 

h  are  in  the  centre-of-mass  system. 

•j* 

The  kinematic  expressions  used  in  the  subroutine  are  given  below. 


Define 


m1m4(e1/et) 


(m1+m2)  (M3+M4) 


MiM3(E!/Et) 


(M1+M2) (M3+M4) 


m2m3 


(M!+M2)  (M3+M4) 


4*  # 

Obtained  from  page  162  of  Nuclear  Data  Tables,  Part  3,  United  States 
Atomic  Energy  Commission  (1960)  . 
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D  = 


M  M 
2U4 


Mi 


(M1+M2)(M3+M4) 


(1  + 


m2 


E? 


where 

Mi  =  mass  of  the  incident  particle 

M2  =  mass  of  the  target 
M3  =  mass  of  the  outgoing  particle 
M4  =  mass  of  the  residual  nucleus 
Ei  =  energy  of  the  incident  particle 
Q  =  Q-value  to  the  state  in  the  residual  nucleus 
and  Ex=  Ei  +  Q 

The  energy  of  the  outgoing  particle  is  then  given  by: 

E3  =  EX  B[cos6lab  +  (D/B  -  sin20lab) 1 7 2 ] 2 

The  positive  sign  is  used  in  the  program  to  reduce  the  computer  running 
time,  since  the  double  values  only  occur  when  the  neutron  energy  is  <  100 
keV.  The  centre-of-mass  angle,  9  ,  of  the  outgoing  particle  is  given  by 

E3  1/2 

cos  6  =  (—  -  B  -  D) /  [2 (AC)  ] 

cm  E^ 

For  the  case  of  n-^He  scattering,  the  centre-of-mass  angle  is  obtained  by 
putting  in  the  above  expression,  the  appropriate  masses  and  setting  Q  equal 
to  zero. 

Al. 2  Description  of  the  Computer  Program 

An  alphabetic  list  of  the  definitions  and  symbols  used  in  the  program 
is  given  in  Table  Al-1.  Tables  Al-2  and  Al-3  give  the  input  and  output 
quantities.  Also  at  the  end  of  the  Appendix,  a  flow  diagram  and  a  Fortran 
listing  are  given. 
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TABLE  Al-1 

An  Alphabetic  List  of  Definitions  and  Symbols 


AMD  mass  of  deuteron 

AMN  mass  of  neutron 

AMP  mass  of  residual  nucleus 

AMT  mass  of  target  nucleus 

AVE(P2)  average  value  of  P2 

AVE(P2C0SPHI)  average  value  of  P2Cos<j> 

CPHI  cos  <J> 

DED  increment  in  deuteron  energy 

DTHE10  increment  in  0 1 q 

E3  energy  of  outgoing  neutron  in  the  finite  geometry  calculation 

EDI  initial  value  of  deuteron  energy 

EDF  final  value  of  deuteron  energy 

EN  energy  of  outgoing  neutron  from  the  nuclear  reaction 

EXN  excitation  energy  of  a  state  in  the  residual  nucleus 

NED  number  of  steps  in  the  increment  of  ed 

NTHE10  number  of  steps  in  the  increment  of  0 1 q 

PPX3  a  running  index  for  increment  in  PX3 

PPY3  a  running  index  for  increment  in  PY3 

PPZ3  a  running  index  for  increment  in  PZ3 

PPX5  a  running  index  for  increment  in  PX5 

PPY5  a  running  index  for  increment  in  PY5 

PPZ5  a  running  index  for  increment  in  PZ5 

PX3  number  of  parts  into  which  the  scattering  volume  is  divided 

along  the  X3~axis 

PY3  number  of  parts  into  which  the  scattering  volume  is  divided 

along  the  Y 3-axis 

PZ3  number  of  parts  into  which  the  scattering  volume  is  divided 

along  the  Z 3-axis 
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Table  Al-1  (continued) 


PX5 

number  of  parts  into  which  the  neutron  detector  is  divided  . 
along  the  X5~axis 

PY5 

number  of  parts  into  which  the  neutron  detector  is  divided 
along  the  Y5-axis 

PZ5 

number  of  parts  into  which  the  neutron  detector  is  divided 
along  the  Z5-axis 

P2X 

P2CJ 

P2X0 

intermediate  value  of  P20  in  the  process  of  summation 

P2XCPH 

P2CCOS(|> 

P2XCP0 

intermediate  value  of  ^2°  coscf>  in  the  process  of  summation 

Q 

Q-value  to  the  ground  state  of  the  residual  nucleus  in  the 
nuclear  reaction 

REACT 

title  of  the  nuclear  reaction 

S10 

distance  between  the  target  and  the  centre  of  the  scattering 
volume 

S20 

distance  between  the  centre  of  the  scattering  volume  and  the 
front  face  of  the  neutron  detector 

THE10I 

initial  angle  of  the  outgoing  neutrons  from  the  nuclear 
reaction  for  the  increment  in  0 ]_ q 

THE10F 

final  angle  of  the  outgoing  neutrons  from  the  nuclear 
reaction  for  the  increment  in  0]_g 

THE1L 

angle  of  outgoing  neutrons  in  the  finite  geometry  calculation 

THE20 

angle  between  the  line  joining  the  target  and  the  centre  of 
the  scattering  volume  and  the  line  joining  the  centres  of 
the  scattering  volume  and  the  neutron  detector. 

THE2C 

scattered  neutron  angle  in  centre-of-mass  system 

THE2L 

scattered  neutron  angle  in  laboratory  system 

THETA 

angle  of  outgoing  neutrons  in  the  laboratory  system 

THETB 

angle  of  outgoing  neutrons  in  the  centre-of-mass  system 

X 

n-4He  scattering  cross  section,  i.e.  0 

XO 

intermediate  value  of  0  in  the  summation  process. 
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TABLE  Al-2f 


Input  Quantities 


Card  No. 

Parameters 

Format 

1 

REACT 

IX,  15H 

2 

AMD,  AMT,  AMN,  AMP 

Q,  EXN 

4F9 . 5 ,  2F7 . 4 

3 

EDI,  EDF,  DED 

3F5. 2 

4 

PZ3 ,  PX3 ,  PY3 

PZ5,  PX5 ,  PY5 

6F6ol 

5 

S10,  S20 ,  THE20 

3F7  o  2 

6 

THE10I ,  THE10F, 

DTHE10 

3F7.2 

i 

All  masses 

are  in  amu,  all  energies 

are  in  MeV,  all  lengths  are 

in  cm  and 

all  angles  are  in  degrees „ 

Refer  to  Table  Al-1  for 

explanation  of  parameters . 
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Table  Al-3.  A  Sample  Output  of  the  Computer  Program 


ANALYSING 

POWER  FOR 

LIQUID  HE4  POLARIMETER 

J)P29 

ED= 

5.00MEV 

EXN= 

0 . 0000MEV 

Q=  0 . 5200MEV 

PZ3: 

=  5.0 

PX3= 

5.0  PY3= 

10.0 

PZ5: 

=  5.0 

PX5= 

9.0  PY5= 

9.0 

S10  = 

=  200.00CM 

S20= 

25.00CM  THE20= 

123.00DEG 

THETA 

THETB  EN(MEV) 

TIME (NS /M) 

AVE(P2) 

AVE  (P2C( 

10.00 

10.47 

5.4889 

30.84 

0.9718 

0.9519 

20.00 

20.93 

5.4655 

30.90 

0.9721 

0.9583 

30.00 

31.36 

5.4274 

31.01 

0.9725 

0.9605 

40.00 

41.75 

5.3761 

31.16 

0.9731 

0.9619 

50.00 

52.08 

5.3135 

31.34 

0.9737 

0.9631 

60.00 

62.36 

5.2419 

31.56 

0.9744 

0.9641 

70.00 

72.56 

5.1637 

31.79 

0.9750 

0.9651 

80.00 

82.65 

5.0818 

32.05 

0.9756 

0.9659 

90.00 

92.72 

4.9986 

32.31 

0.9761 

0.9666 

100.00 

102.68 

4.9167 

32.58 

0.9765 

0.9672 

110.00 

112.56 

4.8387 

32.84 

0.9767 

0.9677 

120.00 

122.36 

4.7665 

33.09 

0.9769 

0.9681 

- 


'•v 

id  e.o 

. 

OO.G^I 
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Flow  diagram  for  the  calculation  of  n-4He  analysing  power 


V 


Read 

THE10I, 

THE10F, 

DTHE10 


V 


« 
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V 


\  Write 

\  React 

ft 

\  Write 

\  ED , EXN , Q 

\  Write 

\  PZ3,PX3,PY3, 
\PZ5  ,PX5,PY5 

¥ 

Write 

S10 , S20 , 
THE20 

II 

XX 


II 


Ill 
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VII 


Increase  Eq 


ED=ED  +  DED 


Yes 


Yes 


END 


FORTRAN  SOURCE  LIST 


SOURCE  STATEMENT 
$IBFTC  MAIN 

C  AVERAGE  ANALYSING  POWER  FOR  LIQUID  HE4  POLARIMETER 

1  FORMAT ( IX, 15H  ) 

2  FORMAT (4F9 . 5 ,2F7 . 4) 

3  FORMAT (3F5. 2) 

4  FORMAT (6F6.1) 

5  FORMAT (3F7. 2) 

6  FORMAT (1H1,50HAVERAGE  ANALYSING  POWER  FOR  LIQUID  HE4  POLARIMETER) 

8  FORMAT ( 15X , 3HED= , F5 . 2 , 3HMEV , 7X , 4HEXN=,  F7 . 4 , 3HMEV , 7X , 2HQ= , F7 . 4 , 
13HMEV) 

9  FORMAT ( / , 15X , 4HPZ3= , F6 . 1 , 6X , 4HPX3= , F6 . 1 , 6X , 4HPY3= , F6 . 1) 

10  FORMAT ( 15X , 4HPZ5= , F6 . 1 , 6X, 4HPX5= , F6 . 1 , 6X , 4HPY5= , F6 . 1) 

11  FORMAT ( 15X, 4HS10= , F7 . 2 , 2HCM, 3X, 4HS20= , F7 . 2 , 2HCM, IX , 6HTHE20= , F7 . 2 , 
13HDEG) 

12  FORMAT (/ , 10X, 5HTHETA, 7X, 5HTHETB , 7X, 7HEN(MEV) , 7X, 10HTIME(NS/M) , 7X, 7 
1HAVE(P2) , 7X, 13HAVE(P2COSPHI) ) 

13  FORMAT ( 10X, F7 . 2 , 5X, F7 . 2 , 5X, F8  =  4 , 6X, F7 , 2 , 10X, F7 . 4 , 7X, F7 „ 4) 

READ  (5,1) 

READ ( 5 , 2 ) AMD , AMT , AMN , AMP , Q , EXN 

READ ( 5 , 3) EDI , EDF , DED 

READ ( 5 , 4 ) PZ3 , PX3 , PY3 , PZ5 , PX5 , PY  5 

READ(5, 5) S10 ,S20,THE20 

READ (5,5) THE10I , THE10F, DTHE10 

20  ED=EDI 

NED= (EDF-EDI) /DED+1 . 5 

NTHE10= (THE10F-THE10I) /DTHE10+1 . 5 

21  DO  57  L= 1 , NED 
WRITE (6, 6) 

WRITE  (6,1) 

WRITE ( 6 , 8) ED , EXN , Q 
WRITE (6 , 9) PZ3 , PX3 , PY3 
WRITE(6,10)PZ5,PX5,PY5 
WRITE ( 6 , 11) S10 , S  20 , THE20 
WRITE (6, 12) 

U=0. 0174533 
THETA=THE10I 

22  DO  55  M=1 , NTHE10 
C  CALCULATION  OF  EN 

CTHETA=COS (THETA*U) 
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23  CALL  ENEG ( AMD , AMT , AMN , AMP , ED , Q , EXN , CTHETA , EN , THETB , TIME) 

C  CALCULATION  OF  THEIL  FOR  OUT-GOING  NEUTRON 

PPY3=1. 0 
X0=0 . 0 
P2X0=0 . 0 
P2XCP0=0 .0 

24  Y3=-5. 08+10 .16*(PPY3-0. 5) /PY3 
PPX3=1. 0 

25  X3=-2.54+5.08*(PPX3-0.5)/PX3 
PPZ3=1 . 0 

26  Z3=-2. 54+5.08  *(PPZ3-0. 5) /PZ3 
PPZ5=1 . 0 

27  IF(Z3**2+X3**2. GT. 6 .45) GO  TO  48 

28  Z5=5.08*(PPZ5-0. 5) /PZ5 
PPX5=1. 0 

29  X5=-4 . 445+8. 89*(PPX5-0 . 5) /PX5 
PPY5=1 . 0 

30  Y5=-4. 445+8. 89*(PPY5-0 . 5) /PY5 

31  IF(X5**2+Y5**2.GT.19.84)GO  TO  42 

32  A1=-SIN(THETA*U) 

B 1=0.0 

Cl=COS (THETA*U) 

G=SQRT(X3**2+Y3**2+(S10+Z3) **2 ) 

A2=X3/G 
B2=Y3/G 
C2=$10+Z3)/G 
CTHEIL=A1*A2+B1*B2+C1*C2 

CALCULATION  OF  THE2L  FOR  SCATTERED  NEUTRON  AND  COS (PHI) IN  LAB. 
SYSTEM 

33  Pl=X5*COS (THE20*U)+(S20+Z5) *SIN(THE20*U) 

P2=Y5 

P3=-X5*SIN(THE20*U)+(S20+Z5) *COS (THE20*U)+S10 
E=SQRT( (P1-X3) **2+(P2-Y3) **2+(P3-S10-Z3) **2) 

A4= (P1-X3) /E 
B4= (P2-Y3) /E 
C4=  (P3-S10-Z3) /E 
CTHE2L=A2*A4+B2*B4+C2*C4 
IF(THETA. EQ . 0. 0) GO  TO  60 
IF (THETA. EQ . 180. 0) GO  TO  60 

34  D=SQRT((B1*C2-B2*C1)**2+(A2*C1-A1*C2)**2+(A1*B2-A2*B1)**2) 

A3= (B1*C2-B2*C1)/D 
B3=(A2*C1-A1*C2)/D 
C3=(A1*B2-A2*B1)/D 

F=SQRT ( (B2*C4-B4*C2) **2+(A4*C2-A2*C4) **2+(A2*B4-A4*B2) **2) 
A5=(B2*C4-B4*C2)/F 
B5=(A4*C2-A2*C4)/F 
C5=(A2*B4-A4*B2)/F 
CPHI=A3 *A5+B3*B 5+C3*C5 
GO  TO  35 
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60  CPHI=1. 0 
C  CALCULATION  OF  E3 

35  CALL  ENEG (AMD , AMT , AMN , AMP , ED , Q , EXN , CTHEIL , E3 , THETB1 , TIME1) 

C  CALCULATION  OF  THE2C  FOR  SCATTERED  NEUTRON  IN  C.M.  SYSTEM 

AMD2=1. 00898611 
AMT2=4. 00387611 
AMN2=1. 00898611 
AMP2=4. 00387611 
Q2=0 . 0 
EXN2=0 . 0 

CALL  ENEG (AMD2 , AMT2 , AMN2 , AMP 2 , E3 , Q2 , EXN 2 , CTHE2L , E32 , THE2C , TIME2) 
C  CALCULATION  OF  PHASE  SHIFTS 

DTM=0 . 0 
DTP=0 . 0 

DZP=U*(180. 0-21. 7684*SQRT (E3) -2 . 63613*E3+. 689548*E3**1. 5) 
IF(E3.LT.8.0)  GO  TO  36 
DTM=U*(-.0725766*E3+.0133187*E3**2) 

DTP=U* ( . 128117*E3- ( . 84296E-02) *E3**2+( . 110985E-02) *E3**3) 

36  IF(E3.GT. 2.0)  GO  TO  38 

37  DOP=U* (684. 93-2312. 49*SQRT(E3)+2493.53*E3-807.582*E3**1. 5) 
DOM=U*(2.45763*E3+3.55942*E3**2-.505754*E3**3) 

TO  TO  39 

38  DOP=U*(. 0648862+180. 513*SQRT(E3) -92. 6748*E3+19.1394*E3**1. 5 
1-1.43251*E3**2) 

DOM=U* (-10. 9749+14 . 8534*E3-. 97544*E3**2+. 0197815*E3**3) 

C  CALCULATION  OF  G  AND  H 

39  CTHE2C=COS (THE2C*U) 

STHE2C=SIN (THE2C*U) 

SDZP=SIN(DZP) 

CDZP=COS (DZP) 

SDOP=SIN(DOP) 

CDOP=COS (DOP) 

S  DOM=  S IN ( DOM) 

CDOM=COS (DOM) 

SDTM=SIN(DTM) 

CDTM=COS (DTM) 

SDTP=SIN (DTP) 

CDTP=COS (DTP) 

HK=(1. 54617E+12) *SQRT(2 .0*AMN*E3) 

HKI=1 , 0/HK 

GRO=HKI*( SDZP*CDZP-SDTP*DCTP*3 . 0/2 . 0-SDTM*DCTM) 

GR1=HKI * ( 2 . 0*SDOP*CDOP+SDOM*CDOM) 

GR2=HKI* ( SDTP*CDTP*9 . 0/2 . 0+3 . 0*SDTM*CDTM) 

GR=GR0+GR1*CTHE2C+GR2*CHTE2C*CTHE2C 

GI0=HKI* (SDZP*SDZP-SDTP*SDTP*3 . 0/ 2 . 0-SDTM*SDTM) 

GI1=HKI*(2 . 0*SDOP*SDOP+SDOM*SDOM) 

GI2=HKI*(SDTP*SDTP*9 . 0 / 2 . 0+3 . 0*SDTM*SDTM) 
GI=GI0+GI1*CTHE2C+GI2*CTHE2C*CTHE2C 
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HR1=HKI * ( SDOP  *CD0P- SDOM*CDOM) 

HR2=HKI * ( SDTP*CDTP-SDTM*CDTM) *3 . 0 / 2 . 0 
HR=HR1*STHE2C+HR2*SIN(2 . 0*THE2C*U) 

HI1=HKI * ( SDOP*SDOP-SDOM*S  EOM) 
HI2=HKI*(SDTP*SDTP-SDTM*SDTM) *3 . 0/ 2 . 0 
HI=HI1*STHE2C+HI2*SIN(2.0*THE2C*U) 
X=GR*GR+GI*GI+HR*HR+HI*HI 
P2X=2 . 0* (GI*HR-GR*HI) 

P  2X=  2 . 0  * ( GI*HR- GR*HI ) 

P2XCPH=P2X*CPHI 

X0=X0+X 

P2X0=P2X0+P2X 

P2XCP0=P2XCP0+P2XCPH 

PPY5=PPY 5+1.0 

IF (PPY5.LT . PY5) GO  TO  30 

PPX5=PPX5+1. 0 

IF(PPX5.LT.PX5) GO  TO  29 

PPZ5=PPZ 5+1.0 

IF(PPZ5.LT.PZ5)  GO  TO  28 

PPZ3=PPZ3+1.0 

IF(PPZ3»LT.PZ3) GO  TO  26 

PPX3=PPX3+1 . 0 

IF(PPX3 . LT.PX3) GO  TO  25 

PPY3=PPY 3+1.0 

IF (PPY3 . LT . PY3) GO  TO  24 

AVP2=P2X0/X0 

AP2CPH=P2XCP0/X0 

WRI TE ( 6 , 1 3 ) THETA , THETB , EN , T IME , AVP  2 , AP  2  CPH 

THETA=THETA+DTHE10 

ED=  ED+DED 

CALL  EXIT 

END 
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FORTRAN  SOURCE  LIST 


SOURCE  STATEMENT 


$IBFTC  ENEG 

C  KINEMATICS  AT  ONE  ANGLE 

SUBROUTINE  ENEG ( AMU1 , AMU2 , AMU3 , AMU4 , El , QO , EXN , C THETA, E3 , THETB , TIME 
1) 

U=0. 0174533 

Q=QO-EXN 

ET=E1+Q 

AMU=  ( AMU1+AMU2)  *  ( AMU3+AMU4) 

A=AMU1*AMU4*E1/ (ET*AMU) 

B=AMU1*AMU3*E1/ (ET*AMU) 

C=AMU2 *AMU3*(  1 . 0+AMUl*Q/  ( AMU2 *ET)  )  /  AMU 
D=AMU2*AMU4*(le 0+AMUl*Q/ (AMU2*ET) ) /AMU 
E3=ET*B*( CTHETA+SQRT  (D/B- (1 . 0-CTHETA*CTHETA) ) ) **2 
CTHETB= (E3/ET-B-D) / ( 2  0 0*SQRT (A*C) ) 

IF(CTHETB «,  LT. (-1.0)) GO  TO  13 
10  THETB=ARCOS(CTHETB) 

GO  TO  12 

13  THETB=3o 14159 
12  THETB=THETB/U 

TIME=72. 24618/SQRT (E3) 

RETURN 

END 


.  .. 


APPENDIX  2 


SOME  SCHEMATIC  DRAWINGS 
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A2-1 


A  schematic  drawing  for  mounting  a  RCA  8575  phototube  to  the 
tail  section  of  the  cryostat 
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A2-2 


A  schematic  drawing  of  the  channel  iron  cart  for  mounting  the 
polarization  system 


TOP  VIEW 


BACK  VIEW 
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